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FOREWORD 
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of  California,  under  USAF  Contract  No.  AF  33(6l6)-6910.  Hie  work  was  administered 
under  the  direction  of  the  Themomechanics  Research  Branch,  Aeronautical  Research 
Laboratory.  The  authors  are  indebted  to  J.  Frisch,  Associate  Professor  of 
Mechanical  Engineering,  for  supervision  of  the  tests  and  analysis  of  the  results 
described  in  Part  III. 

This  report  covers  research  conducted  from  1  February  i960  through  31 
January  1961. 


ABSTRACT 


A  general  small-deflection  theory  governing  the  elastostatic  extension  and 
flexure  of  thin  laminated  anisotropic  shells  and  plates  is  formulated.  The  plate 
or  shell  structure  may  be  composed  of  an  arbitrary  number  of  bonded  layers,  each 
of  which  may  possess  different  thickness,  orientation,  and/or  orthotropic  elastic 
properties.  Donnell- type  equations  for  cylindrical  shells  and  Poisson-Kirchhoff 
plate  equations  are  explicitly  discussed,  along  with  procedures  for  determining 
stresses  in  an  individual  lamina.  Several  methods  of  solution  of  the  system  of 
equations  governing  extension  and  flexure  of  plates  are  discussed  and  illustrated 
with  examples.  Optimization  of  laminate  configuration  is  treated  briefly.  The 
results  of  a  limited  number  of  crack  propagation  tests  of  flat  plate  aluminum 
foil  laminates  in  uniaxial  tension  are  presented. 
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INTRODUCTION 


Analysis  and  prediction  of  the  mechanical  behavior  of  structural  laminates, 
with  a  view  toward  possible  application  in  aircraft  or  missile  structures  is  con¬ 
tingent  upon  obtaining  basic  information  involving:  (l)  determination  of  mechanical 
properties  of  materials,  (2)  development  of  appropriate  stress  analysis  procedures, 
and  (3)  establishment  of  suitable  criteria  for  strength  analysis. 

This  report  describes  a  theory  of  behavior,  for  flat  plate  and  cylindrical 
shell  structures,  appropriate  to  bonded  laminate-type  construction.  The  possibility 
of  orthotropic  elastic  properties  and  arbitrary  orientation  of  elastic  axes  of 
individual  laminas  is  accounted  for.  The  theory  may  be  extended  to  include  aniso¬ 
tropic  layers  with  only  slight  modification.  Preliminary  experimental  results  of 
crack  initiation  and  propagation  studies  for  aluminum  foil  laminates  in  uniaxial 
tension  are  presented. 
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PART  I:  THEORY  OF  THIN  ANISOTROPIC  PIATES  AND  SHELLS 


INTRODUCTION 

A  general  theory  for  the  flexural  and  extensional  behavior  of  thin  laminated 
anisotropic  shells  is  developed  in  Part  I.  The  theory  is  formulated  within  the 
framework  of  the  classical  shell  theory  predicated  on  the  Kirchhoff-Love  hypothesis. 
In  this  theory  the  effect  of  transverse  shear  deformation  and  transverse  normal 
stress  is  neglected.  Love’s  first  approximation  is  employed  in  the  derivation. 

The  procedures  in  developing  these  governing  equations  are  as  follows.  The  stress- 
strain  relations  for  an  individual  lamina  in  generalized  plane  stress  are  estab¬ 
lished.  A  brief  discussion  on  surface  geometry  pertinent  to  the  derivation  is 
given.  Stress -resultants  and  stress -couples  are  formulated  fey  integration  of  the 
components  of  stress  across  the  thickness  of  the  shell.  The  condition  of  equili¬ 
brium  is  then  imposed  vq>on  the  shell.  The  additional  equation  for  the  compatibility 
of  the  reference  surface  is  given.  Boundary  conditions  associated  with  the  boundary 
value  problem  are  discussed.  This  system  of  general  equations  is  then  specialized 
for  cylindrical  shells  and  flat  plates. 

Much  of  the  early  work  in  laminated  shells  was  devoted  to  symmetric  sandwich- 
type  construction  of  isotropic  materials.  A  theory  of  laminated  orthotropic  shells 
was  developed  by  S.  A.  Ambartsumyan  (l),  who  restricted  the  elastic  axes  of  all  the 
laminas  to  run  parallel  to  the  coordinate  axes.  Other  authors  have  studied  pro¬ 
blems  which  are  governed  by  the  same  type  of  differential  equations  as  for  laminated 
plates  and  shells.  Their  work  will  be  cited.  Other  methods  of  solution,  such  as 
pertubation  and  iteration,  are  also  discussed. 


1.  GENERAL  THEORY  OF  THIN  LAMINATED  ANISOTROPIC  SHELLS 

STRESS-STRAIN  RELATIONS  FOR  A  LAMINA  IN  GENERALIZED  PLANE  STRESS 

To  study  laminate  systems,  it  is  first  necessary  to  establish  the  stress- 
strain  relations  for  a  single  lamina.  Consider  an  individual  lamina  whose  middle 
surface  lies  in  the  plane  z  =  0,  with  the  axes  1  and  2  forming  a  right-handed 
orthogonal  coordinate  system.  The  strain-stress  relations  for  a  completely  general 
anisotropic  material  in  matrix  notation  assumes  the  form: 


t  * 
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In  Eq  (l.  1),  (ci,  a q,  ">12)  and  (  £]_.  ^2,  ^  12)  refer  tc  the  average  normal  and 
shearing  stresses  and  strains,  respectively,  over  the  lamina  thickness.  Due  to 
the  symmetry  of  the  compliance  matrix  Sji,  there  are  only  six  independent  constants 
in  generalized  plane  stress  for  an  anisotropic  material.  If  the  material  is  ortho¬ 
tropic  and  the  principal  elastic  axes  coincide  with  the  coordinate  axes,  the 
strain-stress  relations  reduce  to  the  following: 
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The  matrix  of  coefficients  Sqj  in  Eq  (l.2)  is  related  to  the  conventional  elastic 
moduli  and  Poisson's  ratio  in  the  following  manner: 
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It  is  seen  that  the  elastic  properties  of  an  orthotropic  lamina  are  defined  by  four 
independent  constants.  The  factors  2  and  1/2  in  Eq  (1.2)  are  inserted  so  as  to 
make  the  matrix  a  tensor  in  order  that  tensorial  transformations  can  be  carried  out 
in  the  sequel.  Since  most  of  the  further  work  is  devoted  to  orthotropic  materials, 
subsequent  discussion  is  restricted  to  orthotropic  stress-strain  relations.  For 
convenience  of  application,  it  is  frequently  necessary  to  express  Eq  (l.2)  in 
inverse  form,  i.e.,  as  a  stress-strain  relation.  Inversion  of  Eq  (1.2)  gives: 
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STRESS-STRAIN  RELATIONS  FOR  AN  ORTHOTROPIC 
LAMINA  REFERRED  TO  ARBITRARY  AXES 

In  developing  the  theory  associated  with  laminates  composed  of  individual 
laminas  in  which  the  elastic  axes  are  oriented  at  various  angles  relative  to  the 
major  axis  of  the  laminate  proper,  it  is  necessary  to  express  the  stress-strain 
relations  for  an  individual  lamina  referred  to  orthogonal  axes  making  an  arbitrary 
angle  relative  to  the  elastic  axes  of  the  lamina.  Referring  to  Figure  1. 1,  the 
elastic  axes  1,2,  of  the  lamina  are  rotated  through  a  positive  angle  6  relative  to 
the  arbitrary  reference  axes  a,  p. 
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Figure  1. 1 


It  is  desired  to  express  the  stress-strain  equations  with  respect  to  axes  e,  p. 
If  the  desired  equations  take  the  form 


^11  °1 2  2^l6 
C12  ^22  2C26 
°16  C26  2°66 

then  the  problem  is  to  express  the  Cj_j  matrix  in  terms  of  the  Dy  matrix  from 
Eq  (l.4)  and  functions  of  the  angle  9.  This  can  be  accomplished  by  utilizing 
appropriate  transformations  for  stress  and  strain  matrices  as  follows :  It  is 
easily  verified  that 
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STRESS-STRAIN  RELATIONS  FOR  AN  ORTHOTOOPIC 
LAMINA  REFERRED  TO  ARBITRARY  AXES 

In  developing  the  theory  associated  with  laminates  composed  of  individual 
laminas  in  which  the  elastic  axes  are  oriented  at  various  angles  relative  to  the 
major  axis  of  the  laminate  proper,  it  is  necessary  to  express  the  stress-strain 
relations  for  an  individual  lamina  referred  to  orthogonal  axes  making  an  arbitrary 
angle  relative  to  the  elastic  axes  of  the  lamina.  Referring  to  Figure  1.1,  the 
elastic  axes  1,2,  of  the  lamina  are  rotated  through  a  positive  angle  0  relative  to 
the  arbitrary  reference  axes  o.,  p. 
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Figure  1.1 


It  is  desired  to  express  the  stress-strain  equations  with  respect  to  axes  c,  p. 
If  the  desired  equations  take  the  form 


(1.6) 


then  the  problem  is  to  express  the  Cij  matrix  in  terms  of  the  C^j  matrix  from 
Eq  (l.U)  surd  functions  of  the  angle  0.  This  can  be  accomplished  by  utilizing 
appropriate  transformations  for  stress  and  strain  matrices  as  follows :  It  is 
easily  verified  that 
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where  the  matrix  T  is  given  by 
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where  [  T  j  _1  denotes  the  inverse  of  ^  T  j  .  Comparing  Eq  (1.6)  and  (1.10)  gives 
[c]  =  [t]  [  C  )(t]  (1.11) 


Performing  the  matrix  operations  indicated  by  Eq  (l.  ll)  leads  to  the  following 
expressions  for  the  elements  of  the  C  matrix. 
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In  the  preceding  equations,  the  elastic  constants  for  stress-strain  relations 
referred  to  arbitrary  axes  have  been  expressed  in  terns  of  the  four  independent 
orthotropic  constants  and  functions  of  the  angle  ©. 

It  is  convenient  in  application  to  deal  with  the  conventional  elastic  moduli 
Ei,  Eg,  G  and  Poisson's  ratios  1 f j,  2)  2  instead  of  the  elastic  constants  Cij.  This 
can  easily  be  accomplished  by  utilizing  Eq  (1.5),  however,  to  simplify  the  final 
expressions  for  the  general  relationships,  let  the  following  definitions  be 
introduced. 


where  k  and  \  are 
arbitrary  parameters 


(1.13) 


The  four  independent  parameters  which  describe  the  material  properties  are  now 
E,  3/  ,  k,  and  Substituting  Eq  (l. 13 )  into  Eq  (1.5)  and  in  turn  into  Eq  (l. 12) 
gives : 

C11  =  l-'k'a/' 1 f  cos^  ©  +  k  sin**  ©  +  (2k^  +  4x)  sin2  9  cos2  ©1  (l.  14) 
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c22 


c,,  = 

oo 


'12 


ci-  = 

lo 


'2o 


-^r 

1-lzV  [ 

•I^[ 

rib[ 


1:  cos^  9  +  sin^  9  +  (2kV  +  4x)  sin^  9  cos2  9| 


\  (sin^  9  -  cos2  9)2  +  (l  -  21: V  +  k)  sin2  9  cos2  9] 


(l  +  k  -  h\)  sin2  ©  cos2  9  +  1:U  (sin^  9  +  cos^  9)1  (l.  14) 


0 


(kV-  1  +  2k)  sin  9  cos^  9  -  (lcl>-  k  +  2l)  sin^  9  cos  9 
(kV-  1  +  2X)  sin^  9  cos  9  -  (kX^  -  k  +  2x)  sin  9  cos^  9 


SURFACE  GEOMETRY  OF  SHELLS 

Let  a  and  p  be  orthogonal  curvilinear  coordinates  which  describe  the  surface 
whose  cartesian  coordinates  are  given  by  the  equations : 

X  =  X  (a,  p) 

Y  =  Y  (a,  p)  (1.15) 

3  -  2  (a,  p) 

The  line  element  on  the  undeformed  surface  is 


ds' 


A2  da2  +  B2  dp2 


where  A  and  B  are  the  surface  metric  coefficients  defined  by 
A2  =  (X,a)2  +  (Y,a)2  +  (Z,Q)2 
B2  =  (X,&)2  +  (Y,p)2  +(S,p)2 


(1.16) 


(1.17) 


The  comma  in  the  subscript  denotes  partial  differentiation.  The  principal  radii 
of  curvature  Rj_  and  R2  are  related  to  the  metric  coefficients  of  the  surface  by 
the  equations  of  Gauss  and  Codazzi : 


AB 

R1R2 


=  0 


(1.13) 
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(1.19) 


Bx 

R 2ha  ~  R 


Let  the  variables  a,  p,  z  be  the  orthogonal  curvilinear  coordinates  of  the 
space  surrounding  the  surface.  The  line  element  in  this  case  is: 


ds 


A2  (1  ♦  f->‘ 


dcT  +  B"  (1  +  |-)  d£c 

r2 


+  dz 


(1.20) 


ASSUMPTIONS  IN  THE  CLASSICAL  THEORY  OF  SHELIS 

The  following  assumptions  (2)  are  made  in  the  classical  theory  of  shells  to 
simplify  the  mathematical  model. 

In  the  analysis  of  beams,  the  assumption  for  displacements  is  that  plane 
sections  remain  plane  before  and  after  deformation.  This  notion  was  extended  to 
plates  by  Kirchhoff  and  to  shells  by  Love.  This  assumption,  as  interpreted  for 
plates  and  shells,  states  that  normals  to  the  undeformed  surface  remain  normal  to 
the  deformed  surface  and  suffer  no  extension. 


Another  assumption  in  the  classical  shell  theory  is  that  the  thickness  of 
the  shell  is  small  in  comparison  to  the  lateral  dimensions.  As  a  consequence,  the 

quantities  and  are  small  compared  to  unity  and  the  variation  of  the  radii  of 

curvature  over  the  shell  thickness  is  neglected. 


The  component  of  stress  normal  to  the  reference  surface  is  small  in  comparison 
to  the  other  components  of  stress.  This  assumption  stipulates  a  state  of  generaliz¬ 
ed  plane  stress. 


Also  the  strains  and  displacements  of  second  or  higher  order  are  neglected 
in  comparison  to  their  first  order  terms  in  the  classical  theory. 


THE  DISPLACEMENT  VECTOR 

Quantitatively,  the  Kirchhoff -Love  hypothesis  implies  that  the  components 
of  the  displacement  vector  in  orthogonal  curvilinear  coordinates  are: 

w,  u 

u  (a,  p,  z)  =  uQ  (a,  p)  -  z 

v  (a,  p,  z)  *=  vq  (a,  p)  -  (-^  -  g2)  z  (1.21) 

w  (a,  p,  z)  =  wq  (a,  p) 
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In  the  above  equations  uQ,  vQ,  and  vQ  are  the  displacement  components  of  a  point 
on  an  arbitrary  reference  surface. 

STRAIIi-DISPLACStEtlT  3QUATI0IJS 

The  strain-displacement  equations  for  the  linear  theory  of  shells  (also  known 
as  Love's  first  approximation)  are: 


=  e 


Vap  = 


-  sY 

10  ^1 

^20  “  Z  ^2 


(1.22) 


^120  2Z  ^  12 


'Xi,'X2,X 


where  tpo,  *2Q>  *120  are  the  reference  surface  strains  and  A-2>  S*  12 

are  the  chances  of  curvature  and  the  twist  of  the  surface.  The  reference  surface 
strains  are  Given  by  the  same  expressions  for  the  middle  surface  strains  in  the 
theory  of  single-layer  shells.  These  expressions  are  given,  for  example,  by 
Ilovozhilov  (3). 

.  v 

/■  __  X  ,  O  «  .  W 

ei0  =  A  Uo,a  AB  A,p 


Z'  s  —  V"  +  ^  T3  _i_  ^ 

*20  B  o,p  AB  'a  IV 


(1. 23) 


>*  «  -  v  +  -  u 

0 120  A  o,Ci  B  o,p 


AB  ^uo  +  vo 


The  expressions  for  the  changes  of  curvature  which  are  also  given  by  Ilovozhilov 
(3)  are: 


V  _  1  /I  %  1  .  / 1  “on 

A  'A^’a  ■  R/’a  AB  a,b  w,b  "  R.J 


^2  ~  B  w’e  "  R  J'B  +  AB  B,a  ^A  w'a  “  R  ^ 


¥ 

(1.24) 


These  e:g?ressions  ignore  the  extensional  effects  on  the  changes  of  curvi- 
ture.  According  to  Ilovozhilov  ( 3 ) ,  the  effects  are  comparable  to  other  effects 
which  are  neglected  in  this  theory.  The  expressions  which  include  extensional 
effects  are  given,  for  example,  by  Vlasov  (4).  For  certain  types  of  shells  it 
appears  that  for  computational  purposes  the  expressions  which  include  the  exten¬ 
sional  effects  are  more  convenient  to  work  with. 


(1.24) 


AB  (W’C£  ■  A  A,p  w,a 


B  B'a  w,p^ 


u 


,I(1U  _  _£  A  ).l(lv 

R1  'B  0,3  AB  ’3;  R2  'A  o,a 


sV 


SORESS-RESULTANTS  AND  STRESS -COUPLES 


If  the  state  of  generalized  plane  stress  is  assumed  to  exist  in  the  "k-th" 
lamina  of  the  shell,  the  stress-strain  relation  for  this  condition  is  given  by 
Eq  (1.6).  This  equation  rewritten  for  the  "k-th"  lamina  is 


°c(k)  ' 

c  M 
C11 

c  <k> 
C12 

c  (k) 

Cl£ 

°6<k) 

S 

c  <k> 
C12 

c  <k> 

t22 

c  M 
C26 

T  M 

c  (k) 

p  (k) 

p  (k) 

ap 

k 

cl6 

1 

C26 

C66 

i 

(1.25) 


Hie  constants  can  be  found  for  each  lamina  from  Eq  (l.  12)  or  (l.  14)  hy  sub¬ 

stituting  the  value  of  appropriate  to  the  particular  lamina  "k. " 


Stress-resultants  and  stress -couples  can  be  formulated  in  terms  of  the  dis¬ 
placements  by  integrating  Eq  ( 1. 25 )  across  each  lamina  and  sunning  the  resulting 
expressions  over  n  layers. 


V  V 


£  f*  (“a'10’  V°>(1  +  fc> 


k=l 


V] 


dz 


k=l 

Vi 

Hr’ 

n 

I 

t 

(»a<k)-  Vk>>  ^  (1 

k=l 

Vi 

V 

V " 

n 

I 

f* 

J 

Vk>)  1  (1 

k»l 

hk-l 

(1.26) 


(1.27) 


According  to  Love's  first  approximation,  the  quantities  5-  and  are  neglected. 
Therefore  =  N^q  and  . 
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Na 

A11  A12  ai6 

^10 

°11  Dll  hi 

*1 

- 

A12  A22  ^26 

^  20 

- 

DJ  ht  D2? 

*2 

V 

Al6  A2  6  A66 

V120 

hi  D2?  hi 

k 

(1.28) 


M 

a 

'  Dll 

Dl2 

di i 

^10 

D11 

D12 

di6' 

V 

- 

DlS 

°22 

hi 

^20 

- 

°12 

D22 

D26 

*2 

V 

Dl t 

°2^ 

hi 

"^120 

D16 

D26 

D66 

* 

**12, 

(1.29) 


where  the 


■ft 

A±j,  D^j,  and  D.^  are  defined  as 
Aij  *  )  <hk  "  ^-1^ 


°f  j  -  i 

r - 1 

) 

k=l 

5«(k)  ■  ^-i) 

DiJ  *  3 

£ 

k=l 

Eu<k) 

(1-30) 


EQUATIONS  OF  EQUILIBRIUM  AND  COMPATIBILITY 

The  five  equations  of  equilibrium  for  a  shell  element  and  the  compatibility 
equation  for  the  in-plane  strain  components  of  the  reference  surface  constitute 
the  determinative  system  for  this  problem.  These  equations  referred  to  the 
reference  surface  are : 


(BBa>-a  *  <AV'f>  +  "as  %  *  Np  B’a  *  If  1  “  %  ’  0 
<>V-a  -  <««„>,„  *  %  B-a  ’  ”0  %  +  f2  %  *  “  1f>  ’  ° 
<“t>-o  *  <*V'B  '  “  +  +  “  ■  0 

*  (AV'e  *  %  %  -  %  B-o  '  “  %  '  0 


(1.31) 
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(1-31) 


<  *  <V'3 

+  ih£ 

R2  A3l 

aIB 

+  A,  „ 

5  1  10,|3 

€cO  "  6io)  -  7  *i2o-%*J, 


AB  [b  {.A^10,p  +  A,p  ^  6 10  "  ^ 20 ^  '  2  ^120,a  “  B'a  ^12ol  , 


=  0 


(1-32) 


In  the  above  equations  q^,  qp,  and  q2  are  the  external  loads  on  the  shell;  ^ 
and  Qp  are  the  transverse  shearing  forces,  and  all  other  symbols  have  been 
previously  defined. 


BOUNDARY  COnDITIOKS 

Associated  uith  the  system  of  Eq  (l.3l)  and  (1.32)  are  i'our  boundary  condi¬ 
tions  to  be  satisfied  on  the  contour  of  the  shell.  These  boundary  conditions 
arise  from  the  physical  requirements  of  support  along  the  shell  contour.  Expres¬ 
sions  of  these  conditions  may  be  conveniently  obtained  by  variational  methods. 
Such  an  analysis  leads  to  boundary  conditions  which  are  identical  to  those  of 
the  usual  isotropic  shell.  However,  it  should  be  emphasized  that  in  the  present 
case  all  boundary  forces  and  displacements  are  referred  to  the  reference  surface. 
If,  in  the  physical  problem,  the  boundary  forces  and  displacements  are  given  for 
some  other  surface,  an  equivalent  refei-ence  surface  system  must  be  calculated. 

The  boundary  conditions  for  shells  are  summarized  in  the  following  table. 


I 

11 

--.I 

Force  Boundary  Conditions 

Displacement  Boundary  Conditions 

■ 

(N  -  Ti) 
n  n 

(u  -  u  ) 
n  u 

2 

/„  fInt  —  1  ntv 

(Nnt  +  R“  -  *nt  ~  — > 
n  ii 

(ut  ’  “t} 

3 

(i_  w  -  J1  ) 

'A  n  K  1  ' 

n  n 

■ 

^‘r  Mnt,  t  -  V 

_ a _ _ _ 

(w  -  w) 

(i.  33) 
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where 

Nn  -  normal  force 

Nt  -  in-surface  shearing  force 

Qjj  -  transverse  shearing  force 

Mn  -  bending  moment 

Mot  -  twisting  moment 

Vn  -  effective  shear  given  by 

V  =  Q+v-M,. 
n  ti  An  nt,  t 

S  -  nonaal  displacement 
u^  -  tangential  displacement 
w  -  deflection  normal  to  surface 

-  angle  of  rotation  in  the  normal  direction  given  by 


An  -  metric  coefficient  associated  with  the  contour. 

Bars  over  the  quantities  signify  that  they  are  the  prescribed  values  at  the 
boundaries . 

In  each  of  the  four  conditions,  either  the  barred  quantity  in  column  I  or  II 
must  be  prescribed.  For  that  force  or  displacement  which  is  prescribed,  the 
bracketed  quantity  must  vanish.  Uxis  means  that  the  unbarred  quantity  (that  force 
or  displacement  from  the  interior  of  the  shell)  must  take  on  the  prescribed  value 
on  the  boundary. 

Some  examples  of  boundary  conditions  are : 

Free  edge : 

M  . 

N  =  0  ,  I!  +  -25  c  o  , 
n  nt  R  ’ 

n 

Hinged  edge  with  fixed  support: 

M  =  0  ,  u  -  0  ,  u  =  0  ,  w  »  0 

XX  U  u 

Hinged  edge  with  support  free  to  move  in  the  normal  direction: 

M  =  0  ,  w  =  0  ,  u  =  0  ,  H=0 

xi  t  n 
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Clamped  edge: 


% 


0  , 


VL 


0  , 


w  ■  0  , 


0 


The  theory  developed  need  not  be  restricted  to  shells  comprised  of  ortho¬ 
tropic  laminae.  Because  of  possible  rotation  of  the  orthotropic  elastic  axes, 
the  equations  already  involve  six  Cit»,  which,  however,  are  not  independent  for 
a  shell  of  orthotropic_laminas .  In  the  most  general  case,  there  can  exist  six 
independent  constants  C^j.  It  is  then  only  necessary  to  regard  these  six  con¬ 
stants  as  completely  independent  for  a  theory  which  involves  a  completely  aniso¬ 
tropic  material. 
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2.  DONNELL-TYPE  EQUATIONS  FOR  LAMINATED  ANISOTROPIC  CYLINDRICAL  SHELLS 

In  this  section,  the  general  equations  from  the  preceding  section  are 
specialized  for  a  circular  cylindrical  shell  using  the  Donnell  approximations, 
which  neglect  certain  tenns  in  the  changes  of  curvature  and  in  the  equilibrium 
equations.  The  justification  for  neglecting  these  terns  can  be  argued  by  con¬ 
sideration  of  the  configuration  of  the  shell  and  the  geometry  of  defomation. 

The  Donnell  Approximations  simplify  the  mathematical  model  considerably  and  the 
results  obtained  from  this  theory,  as  discussed  in  connection  with  isotropic 
cylindrical  shells  by  Kempner  (5),  agree  closely  to  the  results  obtained  from  a 
more  exact  theory  in  most  cases  of  technical  interest. 

DERIVATION  OF  THE  GOVERNING  DIFFERENTIAL  EQUATIONS  AND  BOUNDARY  CONDITIONS 

The  cartesian  coordinates  which  describe  a  circular  cylinder  of  radius  a 

are : 

X  =  a  cos 

Y  =  a  sin«f  (2.1) 

Z  =  x 

Eq  (2.1)  has  the  form  of  Eq  (l.  15)  with  the  variables  a  =  x  and  p  *  .  Die 

metric  coefficients  can  be  computed  frcm  Eq  (l. 17).  Die  square  of  the  line 
element  in  this  system  of  curvilinear  coordinates  is  then 
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Stress -resultants  and  stress-couples  are  formulated  as  described  in  the 
previous  section.  These  expressions  for  a  cylindrical  shell  are : 


N 

X 

'*n 

A1 2 

ai6 

u 

o,x 

* 

'  - 

A12 

^22 

*26 

1  w 

-  V  +  - 

a  o,f  a 

N  , 

Al£ 

► 

^6 

^66 

+x 

o 

> 

M 

X 

'D1 1 

Dll 

Dl2 

uo,x 

Mp 

a 

Di! 

^22 

D2? 

iv  +2 

a  o/f  a 

Di£ 

°2? 

D<£, 

v  +  ■iu 
o,x  a  o ,tf 

D11  Dll  D1 % 

'  W>xx 

Dll  D2l  D2? 

a 

Dl2  D2?  °6b, 

2 

— w,  . 
a  X^ 

(2.5) 


D11 

D12 

Di6 

w, 

XX 

°12 

D22 

D26 

1 

a 

di6 

D26 

D66 

2 

-w, 

a  ’x^ 

(2.6) 


The  five  equilibrium  equations,  which  are  also  modified  by  a  Donnell  approxi¬ 
mation,  supplemented  by  the  compatibility  condition  for  the  deformed  reference 
surface  of  the  shell  constitute  the  determinative  system  of  equations  for  the 
problem.  These  equations  are: 


Nx,x  +  5  %,* 


Sc,x  +  a  -  S  + 


q  ■  0 
z 


M  +  -  M  ..  -  Q  «  0 

x,  x  a  xif^p  x 


(2.7) 


M  +  -  _  -  Ss  =  0 

xy,x  a  *,*  * 

a  w,xx  +  ^20, xx  +  ~2  '  a  ^120,xtf  "  0 

cL 


In  the  equilibrium  Eq  (2.7),  the  only  external  load  which  is  considered  is  the 
normal  load  qz-  0*  and  appearing  in  the  above  equations  can  be  eliminated  by 


appropriate  differentiation  and  subsequent  addition  of  the  fourth  and  fifth 
equilibrium  equations  using  the  third  equilibrium  equation  as  an  identity. 


M  +  2-  M 
x,  xx  a  x«f,xtf 


M 


0 


(2.9) 


It  is  convenient  to  treat  the  transverse  deflection  w  and  the  8 tress - resultants 
Nx*  M  priniary  dependent  variables  of  the  problem.  Equations  (2.5) 
and  (d.6)  rewritten  in  abbreviated  matrix  notation  are: 

N«A€0-D*%  (2.5a) 


M  -  D*€  -  D7C 

o 


(2.6a) 


Inversion  of  Eq  (2.5a)  gives: 

€  -  BN  +  BD* 

'"o 


(2.10) 


where  B  is  the  inverse  of  A. 


B11 

B12 

B12 

B22 

B26 

Bl6 

B26 

B66 

(2.11) 


The  components  of  the  symmetric  B  matrix  are : 

A12A26  ‘  ^^22 
A12Al6  ' 

where  ^  ls  "the  determinant 


An 

A12 

ai6 

^A  - 

A12 

Ag2 

*26 

ai6 

A26 

A66 

"  A26 


Al6A26  "  A12A66 


aua66  "  ai6 
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Substituting  Eq  (2.10)  into  Eq  (2.6a)  gives: 


dl6  ~  DllDl6Bll  +  D12Dl6B12  +  Di^Bl6  +  D11D26B12  +  D12D26B22 

+  Dl6D26B26  +  DUD66Bij6  +  D12D66B26  +  Dl6D66B66  “  Dl6 

**22  D12  B11  +  2DI2D22B12  +  2D12D26Bl£  +  °22  B22  +  2D22D26B26  +  D26  B66  D22 

d26  ”  Dl2°l6Bll  +  Dl6D22B12  +  Dl6D26Bl6  +  D12D26B12  +  °22D26B22 

+  D262b26  +  D12D66Bl6  +  D22D66B26  +  D26D66B66  “  D26 

d66  "  D!62bii  +  +  2Dl6D66Bl6  +  D262b22  +  2D26D66B26  +  D6 6B66  "  D66 


way  of  recapitulation,  Eq.  (2.10)  and  (2.12)  are  restated  in  terms  of  the 
components  of  the  property  matrices  defined  previously. 

«10  ■  “o,x  ■  BlA  +  BlA*  +  Bl6"x»  ♦  ‘■ii1'.**  +  h  '"f»  +  2bk  s 

€20  ‘  i  To,„  -  I  -  V*  +  B22"r  +  B26M*f  +  b2l“'«  +  »k 


+  2bjL;  -  w, 

26  a  x«f 


(2.10a) 


*120  ■  «„,»  ■  BlA  +  %6V  Wif*  b6l“-xx  +  b62  p  »•*» 

♦a*... 

a  'x«f 

Mx  *  biA  *  V*  -  bi6V  *  duw,«x  +  ~jr  v"ft  *  -T5  “•** 

d22  2dp6 

Mf  “  b21*x  +  b22Nf  +  b26Nxf  +  d12W,xx  +  +  V’x<f 

d?6  24^- 

Mx^  "  b6lNx  +  b62«f  +  b66Nxf  +  dl6w,xx  +  JT  +  “T6  w'x* 


(2.12a) 


Introduction  of  the  Airy  stress  function  U  defined  in  cylindrical  coordinates 
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1_  u 
2 

a 


lu'*f 


(2.15) 


which  identically  satisfie  3  the  in-surface  equilibrium  equations,  the  first 
two  equations  of  Eq  (2.7),  reduces  the  governing  system  to  two  primary  dependent 
variables.  Substitution  of  Eq  (2.10)  into  Eq  (2.8)  using  this  definition  of  the 
Airy  stress  function  gives : 


■“22  u,xxxx  a  u,xxxtf  q2  ’xxff 


»3  'xifff 


Bn  (2b,  “ 

♦  -r 

(bu  ♦  ^  -  2^)  .  (»61  -  v 

- y - — y —  •**** 


(2.16) 


+ -|iw, 


-  -  w,  «  0 
a  xx 

Substitution  of  Eq  (2.12)  into  Eq  (2.9)  again  using  Eq  (2.15)  gives: 

b  u,  +(a,& :  V1  o,  +  (Vt  „ 

12  xxxx  a  xxx  if  2  xxifif 


.  (ab6l  -  b26)  .. 
+  - n -  u, 


+  1  ’i  *  ■  Uj iMUHMM  d.  -  V.  +  V. 

**■  ,,fWf  11  xxxx  a  xxxtp 


.  2<d12  +  2d66>  _r  .  "d26  lr  .  d22  lr 
+ -  -  ’xx^f  +-XW’xW 


(2.17) 


-  U,  »-q 
a  xx  Hz 


Equations  (2.16)  said  (2.17)  constitute  the  generalization  of  the  Donnell 
equations  of  cylindrical  shells  to  include  the  effect  of  anisotropy  of  the 
individual  Isminas.  Hie  presence  of  both  dependent  variables  indicates  coupling 
between  membrane  and  bending  effects. 

The  boundary  conditions  are  given  by  Eq  (1.33).  Hiey  are  summarized  below 
to  complete  the  formulation  of  the  cylindrical  shell  problem. 
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Force  Boundary  Conditions 


Displacement  Boundary  Cond. 


(N  -  N  ) 
v  n  n7 


M  .  M  . 

(N  +  ■  -  N _ —) 

v  nt  R  nt  R  ' 

n  n 


(M  -  M  ) 
n  n 


(u  -  u  ) 
n  n' 


<ut  -  V 


(v  -  w) 


/ 


DETERMINATION  OF  STRESSES  IN  AN  INDIVIDUAL  LAMINA  AND 
INTER-LAMINAR  SHEAR  STRESSES 

The  stresses  in  the  "k-th"  lamina  of  a  shell  are  given  by  Eq  (1.6).  They 
can  he  written  in  matrix  form  using  the  primary  dependent  variables  U  and  w  by- 
substituting  Eq  (2.10)  into  Eq  (l.22)  and  in  turn  into  Eq  (1.6).  Die  expression 
for  stresses  in  matrix  notation  then  becomes : 


where 


and 
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p(k)B  4p(k)R  +c(k)B  7?(k)B  ^k)B  ^(k) 
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It  should  be  noted  that  Tv  '  and  t'  '  are  unsyrametric.  The  stresses  written 
explicitly  are : 
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Although  transverse  shear  deformation  is  neglected  in  the  theory  discussed, 
transverse  shear  resultants  can  be  determined  from  the  equilibrium  equations, 
i.e.  the  last  two  equations  of  Eq  (2.7),  and  are  defined  as: 


n  \ 

Q*  -  V  f  T^zk)  dz  (2-22) 
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In  terms  of  the  primary  dependent  variables  U  and  w,  these  transverse  shear 
resultants  are : 
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The  inter- laminar  shear  stresses  t  and  t . are  formed  by  consideration 

xz  *f  z 

of  equilibrium  in  the  x  and  «f  directions.  Summing  the  forces  in  the  x-direction 
up  to  the  (k-th)  layer  gives : 


k  h 
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Here  fxzo  and  t  ^  zo  are  constants  to  be  adjusted  by  conditions  at  the  top  and 
bottom  shell  surfaces.  Using  Eq  (2.19a),  these  expressions  become 
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3-  EQUATIONS  FOR  LAMINATED  ANISOTROPIC  PLATES 


An  Important  area  of  application  is  that  of  laminated  anisotropic  plates. 
In  this  section,  the  general  equations  of  Section  1  are  specialized  for  rectan¬ 
gular  plates. 

DERIVATION  OF  THE  GOVERNING  DIFFERENTIAL  EQUATIONS  AND  BOUNDARY  CONDITIONS 


T-  x 


* 


Figure  3-1 


A  plate  is  defined  as  a  shell  with  no  curvature;  therefore,  Ri  “  Rg  “  ® 
and  the  line  element  is  given  try: 

ds2  =  dx2  +  dy2  (3-1) 

where  the  change  of  variables  of  a  ■  x  and  p  *  y  from  Eq  ( 1. 16)  has  been  made. 

The  reference  surface  strains  given  by  Eq  (l.23)  reduce  to  the  following 
for  a  plate: 
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(3-2) 


The  expressions  for  the  changes  of  curvature  are  now: 


(3-3) 
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The  3 tress -resultants  and  stress -couples,  defined  by  Eq  (1.26)  and  Eq  (1.27), 
for  rectangular  plates  become: 
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are  defined  by  Eq  ( 1. 30 ). 

(3-M 


(3-5) 


The  three  equations  of  equilibrium  for  the  plate  element  supplemented  by 
the  compatibility  equation  for  the  in-plane  strain  components  at  the  reference 
surface  of  the  plate  constitute  the  determinative  system  of  equations.  These 
equations,  as  reduced  from  Eq  (l.  31)  and  Eq  (1.32),  are 
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It  is  convenient  to  regard  the  transverse  displacement  v  and  the  stress- 
resultants  Nx,  Ny,  and  as  the  primary  dependent  variables  for  this  system 
of  equations.  The  moment  and  reference  surface  strains  can  be  expressed  in  terms 
of  these  variables  in  the  same  manner  as  for  cylindrical  shells. 
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where  B,  b,  b  ,  and  d  are  defined  hy  Eqs  (2.11),  (2. 13),  and  (2.14).  Eq  (3.9) 
and  Eq  (3*10)  written  explicitly  are: 

£10  -  “o,x  "  BlA  +  Vy  +  Bl6»xy  +  h^xx  +  4v'yy  +  ^’xy 

*20  ‘  'ey  -  BlA  +  W  +  Vxy  +  b2lv'xx  +  V'-yy  *  ^-xy  '3'9“> 

*120  *  uo,y  *  vo,x  -  BlA  *  B26Hy  +  B66Hxy  +  b6lv'xx  +  b&>v'yy  *  ^''xy 


M  =  b.-.N  +  b10N  +  b.^N  +  d,  ..w,  +  a„0v,  +  2d^w, 

x  11  x  12  y  lb  xy  11  'xx  12  yy  lb  'xy 

M  =  b01N  +  b__N  +  b  ,N  +  d_0W,  +  cU_w,  +  2d^-w, 

y  21  x  22  y  26  xy  12.  xx  ^22  yy  ^6  xy 
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M  =  b,,N  +  V_K  +  b/-/-N  +  d.y-w,  +  cLxW,  +  2A,cv, 

xy  61  x  62  y  66  xy  16  'xx  i26  'yy  66  'xy 


Introduction  of  the  Airy  stress  function  U  defined  in  rectangular  cartesian 
coordinates  by: 
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which  identically  satisfies  Eq  (3*6)  reduces  the  determinative  equations  to  two 
primary  dependent  variables.  Substituting  Eq  (3-9)  into  Eq  (3*8)  and  using  Eq 
(3-11)  gives: 
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+  ^2b6l  ~  b26^W'xyyy  +  b21W'yyyy  *  0 
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Substituting  Eq  (3.10)  into  Eq  (3-7)  again  using  Eq  (3. 11)  gives 


b10U,  +  (2b.0  -  b.^)U,  +  (b,.,  +  b00  -  2b^)U, 

12  ’  xxxx  v  to2  16'  xxxy  11  22  06'  ’xxyy 
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ol  26  xyyy  21  yyyy  11  xxxx  16  'xxxy 

+  2(d10  +  2d/---)w,  +  4d0,-w,  +  d00w,  =  -  q 

'  12  60'  ’xxyy  2 o  xyyy  22  ’yyyy  ^ 


Equations  (3-12)  and  (3-13)  are  the  generalizations  of  the  Poisson-Kirchhoff 
equations  to  take  into  account  the  orthotropic  properties  of  the  laninas  composing 
the  plate.  The  presence  of  both  dependent  variables  in  each  of  the  two  equations 
indicates  coupling  between  the  in-plane  forces  and  the  deflection  surface. 

The  boundary  conditions,  given  by  Eq  (1.33);  are  summarized  below  to  complete 
the  formulation  of  the  plate  problem. 
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normal  force 
in-p lane  shearing  force 
transverse  shearing  force 
bending  moment 
twisting  moment 
effective  shear  given  by 

V  =  Q  +  M  .  . 
n  tl  nt,  t 
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%  -  normal  displacement 

ut  _  tangential  displacement 

v  -  transverse  deflection 

w,n  -  slope  in  the  normal  direction. 

Bars  over  the  quantities  signify  that  they  are  the  prescribed  values  at  the 
boundaries . 

In  each  of  the  four  conditions,  either  the  barred  quantity  in  column  I  or 
II  must  be  prescribed.  For  that  force  or  displacement  which  is  prescribed,  the 
bracketed  quantity  must  vanish.  Ibis  means  that  the  unbarred  quantity  (that  force 
or  displacement  frcm  the  interior  of  the  plate)  must  take  on  the  prescribed  value 
on  the  boundary. 

The  boundary  conditions  corresponding  to  the  physical  restraint  may  be 
thought  of  in  terns  of  those  of  the  usual  conditions  for  a  bending  problem 
(conditions  3  and  4)  and  an  in-plane  problem  (conditions  1  and  2).  Certain 
plate  problems  may  be  decomposed  into  these  two  separate  problems  and  treated 
individually.  However,  in  general  the  problem  must  be  treated  as  a  combined 
bending  and  in-plane  problem  since  U  and  w  are  coupled  in  the  governing  differ¬ 
ential  equations  and  in  the  boundary  conditions. 


Some  examples  of  boundary  conditions  for  a  rectangular  plate  are: 


Clamped  Edge  -  In  a  plate  with  a  built-in  edge  both  displacements  Uq  and 
vQ,  the  deflection  w,  and  the  slope  w,n  are  zero.  This  condition  along  the  x 
axis  is 
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(3-15) 


uQ  and  v  may  be  expressed  in  terns 
integration  or  Eq  (3* 9a) 
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of  U  and  w  in  a  specific  problem  by  the 

b,  ,w,  +  b  w,  +  2b£lw,  )  dx  +  f(y) 

11  'xx  21  yy  61  acy  1 

b,0w,  +  b00w,  +  2b_^w,  )  dy  +  f0(x) 

12  'xx  22  'yy  26  xy  2' 


The  arbitrary  functions  fi(y)  and  fg(x)  may  be  evaluated  by  differentiating 
\Xq  with  respect  to  y  and  vQ  with  respect  to  x  and  substituting  into  the  third 
Eq  (3-9a). 
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Free  Edge  -  Along  a  free  edge  of  a  plate,  the  in-plane  and  transverse  forces 
and  the  bending  moment  vanish.  The  displacements  and  slope,  hovever,  are  un¬ 
specified.  Along  the  x-axis,  this  condition  is 
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My  and  Vy  may  be  expressed  in  terms  of  U  and  v  by  Eq  (3. 10a)  and  Eq  (3-21). 
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Simply  Supported  Edge  -  A  simply  supported  edge  permits  rotation  but  no 
deflection.  Since  it  is  free  to  rotate,  the  bending  moment  is  zero.  For  the 
in-plane  boundary  conditions,  either  a  fixed  (pinned)  or  a  free  (roller)  condi¬ 
tion  may  be  prescribed  depending  on  nature  of  the  restraint.  For  a  roller 
support  which  produces  no  in-plane  forces  along  the  x-axis,  the  boundaiy  condi¬ 
tions  are 
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My  in  terms  of  U  and  w  is  given  by  Eq  (3.16a). 


(3-17) 
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DSTSRl-IIKATIOX  CF  STRESSES  IE  AN  INDIVIDUAL  LAMINA  AED 
INTER- LAMINAR  SHEAR  STRESSES 

Tne  general  expressions  for  the  stresses  in  the  ”k-th"  lamina  are  given  in 
Eq  (l.6).  For  the  plate,  this  equation  in  terms  of  the  primary  dependent  vari¬ 
ables  U  and  w  becomes : 
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where  T^k^,  t^k^,  and  C^k^  are  previously  defined  by  Eqs  (2.20),  (2.21),  and 
(1.12),  respectively.  Equation  (3.16)  written  explicitly  are: 
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Although  transverse  shear  deformation  is  neglected  in  this  theory,  transverse 
shear  resultants  can  be  determined  from  equilibrium  considerations.  The  result¬ 
ants  are  defined  as : 
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xz 


dz 
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32 


These  stress  resultants  are  determined,  from  the  moment  equilibrium  equations 
i.e.,  the  last  two  expressions  of  Eq  (l.3l): 
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x,x  xy,y 
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xy,x  y,y 


Combining  these  equations  with  Eq  (3*10)  gives: 
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Figure  3-2 


(  k)  (k ) 

The  inter- laminar  shear  stresses  "Xz  '  and  Tyz  '  are  found  by  consideration 
of  equilibrium  in  the  x  and  y  directions.  Figure  3*2  shows  only  the  change  of 
the  stresses  in  the  x  ar.d  y  directions.  Summing  the  forces  in  the  x-direction 
up  to  the  k-th  layer  gives 

k  h. 
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xz  L*  J  x,x  xy,  y 
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here  *J;zo  and  fyzo  are  constants  to  be  adjusted  by  conditions  at  the  top  and 
bottom  plate  surfaces.  Using  Eq  (3.18a),  these  expressions  become 
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4.  METHODS  OF  SOLUTION 


Several  methods  of  obtaining  solutions  to  the  plate  Eqs  (3.12)  and  (3*13) 
are  discussed  in  this  section.  First,  a  method  which  involves  the  combination 
of  two  fourth  order  equations  into  one  eighth  order  is  indicated.  Then,  iterative 
and  perturbation  schemes  are  discussed  as  other  possible  methods  of  attack.  For 
special  orientations  of  the  individual  laminas,  the  differential  equations  for 
plates  and  shells  simplify,  and  an  extensive  amount  of  information  for  solving 
these  simplified  equations  is  available. 

REDUCTION  OF  THE  COUPLED  EQUATION'S  TO  ONE  EIGHTH  ORDER  EQUATION 

A  possible  method  of  solution  of  the  coupled  system  of  equations  is  to 
combine  the  two  equations  into  one  eighth  order  equation  using  the  technique 
which  Vlasov  used  to  solve  shallow  shell  equations  (7).  For  convenience,  Eqs 
(3-12)  and  (3-13)  are  rewritten  in  operator  notation: 
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4^  (U)  +  L^}  (w)  =  -  q 
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(k.Z) 


(ll)  (  li.  )  (  ll ) 

where  operators  ,  L^  ,  L'  are  the  following  fourth  order  operators: 
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If  there  exists  a  function 


such  that 
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v  =  ii4)  (4>) 
u  =  -  (♦). 

then  Eg  (4.1)  is  satisfied  identically  and  Eq.  (4.2)  becomes: 

43,<*>-  {[4*>]2-4*>4mH  ■< 

where  48)  is  the  following  eighth  order  operator: 
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The  solution  of  Eq  (4.5)  nay  now  be  separated  into  a  homogeneous  and  a 
particular  solution.  The  particular  solution  should  satisfy  the  loading  condi¬ 
tion  but  may  not  necessarily  satisfy  the  boundary  conditions.  This  solution  may 
be  determined  by  the  use  of  trigonometric  series.  The  homogeneous  solution  must 
be  such  that  the  sum  of  the  particular  and  homogeneous  solutions  satisfies  the 
boundary  conditions.  The  latter  condition  may  be  satisfied  in  the  same  manner 
as  that  used  by  M.  Suchar  (6)  to  determine  influence  surfaces  for  anisotropic 
plates,  or  by  other  means. 


ITERATIVE  SOLUTION  OF  THE  COUPLED  EQUATIONS 

Since  Eqs  (3*12)  and  (3*13)  are  based  on  small  deflection  theory,  the 
coupling  between  the  in-plane  stress  function,  U,  and  the  lateral  deflection, 
w,  generally  will  have  a  minor  effect  on  the  solutions.  Hence,  as  a  first 
approximation,  the  system  can  be  treated  an  an  uncoupled  system  by  emitting  from 
the  equations  the  terms  which  contain  the  bji  coefficients.  Then,  by  substituting 
the  first  approximation  into  the  terms  containing  the  b^  coefficients,  a  second 
approximation  may  be  obtained.  The  process  is  then  repeated  to  obtain  higher 
order  approximations. 


Adopting  the  operator  notation  used  previously,  Eqs  (3-12)  and  (3*13)  may 
be  written  as : 


(Vi} 
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(wn) 
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where  n  represents  the  order  of  the  approximation,  wQ  and  U0  being  equal  to  zero. 
It  can  be  seen  that  this  system  is  effectively  uncoupled  since  the  higher  order 
approximations  of  one  variable  are  coupled  to  the  next  lower  order  approximation 
of  the  other  variable.  Hence,  the  methods  of  solution  applicable  to  the  uncoupled 
equations  given  in  the  literature  and  summarized  below  may  be  applied  to  this 
system  of  equations. 

Since  coupling  also  exists  in  the  expressions  for  mcment  and  for  in-plane 
strain,  Eqs  (3*9a)  and  (3.10a),  which  are  used  in  stating  the  boundary  conditions, 
these  expressions  must  also  be  rewritten  in  iterative  form: 
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F3RTUR3ATICN  .SOLUTION  OF  TEE  CO  UP  LSD  EQUATIONS 

For  many  materials  of  practical  importance,  the  individual  laminas  will  be 
only  slightly  orthctropic.  In  such  cases,  the  operators  in  Eqs  (3*12)  and  (3-13) 
will  be  only  slightly  different  from  those  in  the  equations  for  isotropic  plates. 
This  suggests  that  such  problems  may  be  treated  by  a  perturbation  method  of  solu¬ 
tion  (9.  10);  that  is,  the  solution  for  a  substitute  isotropic  material  may  be 
iised  as  a  first  approximation,  and  then  successive  corrections  may  be  obtained 
to  account  for  the  crthotropy. 

To  apply  perturbation  theory,  Eqs  (3-12)  and  (3*13)  are  written  in  a  foim 
such  that  they  may  be  seen  to  be  the  equations  for  an  isotropic  plate  modified 
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where  B  and  D  are  the  maxi  mums  of  %22>  or  ^12  +  ^66>  ^  Dll>  ^22>  or 

Dl2  +  2Q66,  respectively.  This  latter  condition  guarantees  that  the  parameters 
kij,  Jt i j ,  m-y,  nij,  and  n-j?  are  always  le3S  than  1.0  in  magnitude.  Four  con¬ 
stants  have  been  used  to  define  the  substitute  isotropic  material:  B,  2/ 3,  D, 
andVj)  .  This  has  been  done  for  convenience,  since  most  problems  will  involve 
either  an  in-plane  stress  problem  or  a  bending  problem.  The  nunber  of  parameters 
has  also  been  selected  for  convenience  in  the  application  of  the  method. 

When  Eqs  (3-12)  and  (3-13)  are  expressed  in  terms  of  these  parameters  they 
become 
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It  may  be  seen  that  if  the  parameters  kjj,  Xj*,  mij>  nij>  and  n^*  are  so  small 
that  the  terns  containing  then  may  be  neglected,  the  equations  for  isotropic 
plates  are  recovered. 

Since  U  and  v  may  be  considered  to  be  continuous  functions  of  the  parameters 
kij,  jGij,  mi(<,  ni(j,  and  n^,  they  may  be  expanded  in  series  in  teims  of  powers 
of  the  parameters.  However,  from  a  computational  standpoint,  this  method  of 
solution  is  of  practical  value  only  when  the  terns  containing  squares,  products, 
and  higher  powers  of  the  parameters  are  negligible.  Therefore,  the  power  series 
expansions  will  be  used  in  the  following  form: 


U  -  U  +  k,  .  U,  .  .  +  0  . 
o  ij  kij  *ij 

w  =  w  +  k .  .  w.  .  .  +  I.  . 
o  ij  kij  *ij 


U 

w 


+  mij 

Hnij 

+  n.  . 
ij 

UniJ 

U  * 
nij 

+  in  , 

w 

+  n  . 

V 

+  n  * 

w  .  * 

fid 

id 

mij 

id 

nij 

id 

nid 
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The  substitution  of  these  expansions  for  U  and  w  into  Eqs  (4.12)  and  (4.13) 
results  in  the  following  sets  of  equations  when  the  coefficients  of  like  powers 
of  the  parameters  are  equated: 
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Equations  (4.15)  are  the  differential  equations  of  the  first  approximation. 
Equations  (4.16),  (4.17),  (4. 18),  and  (4.19)  are  the  differential  equations  of 
the  corrective  terms  forming  the  second  approximation.  It  should  be  pointed 
out  that  the  equations  involving  w,t  w^jj,  UQi^,  and  Umi .  have  been  emitted 
from  this  summary  since  these  equations  are  all  or  the  formJ  ^(w,  U)  =  0,  and 
their  solutions  are  identically  zero. 

Since  the  elastic  constants  appear  in  the  expressions  for  the  bending 
moments  and  in-plane  strains,  the  boundary  conditions  for  the  above  sets  of 
equations  will  also  be  perturbed.  1/hen  the  bending  moments  are  expressed  in 
terms  of  the  perturbation  parameters,  they  become: 
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These  moments  must  now  be  expanded  in  the  same  manner  as  for  U  and  w,  i.e. 

x  xo  ij,  xklj  Aij  x  JUj  lj  xmij  ij  xnij  ij  xnij 

Substitution  of  the  expansions  for  M,  U,  and  w  into  Eq  (4.20)  leads  to  the  follow¬ 
ing  sets  of  equations : 
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Similar  expressions  may  be  written  for  the  remaining  sets  of  equations. 
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SPECIALIZATION  OF  THE  FLATS  EQUATIONS  FOR  PARTICULAR  ORIENTATIONS 

The  two  governing  Eqs  (3-12)  and  (3-13)  can  be  specialized  for  certain 
orientations  of  the  laminas  composing  the  laminated  plate.  In  the  following, 
the  discussion  is  restricted  to  plates  composed  of  laminas  of  identical  thick¬ 
ness.  Tiro  general  cases  will  be  discussed:  pairwise  orientation  with  symmetry 
about  the  middle  surface  and  spiral  orientation.  Methods  of  solution  associated 
with  each  of  the  cases  are  indicated,  and  wherever  the  equations  have  already 
been  solved,  the  literature  is  cite.. 

PAIRWISE  ORIENTATION  WITH  SYMMETRY  ABOUT  THE  MIDDLE  SURFACE 

This  classification  refers  to  a  laminated  plate  made  from  pairs  of  the  same 
orthotropic  laminas.  I3ie  laminas  are  situated  such  that  one  lamina  of  a  pair  is 
at  the  same  distance  above  as  the  other  is  below  the  middle  surface  (which  is 
also  taken  as  reference  surface).  The  principal  axes  of  both  laminas  of  the  pair 
also  run  in  the  same  direction. 

This  orientation  uncouples  the  in-plane  and  bending  effects  since  D^j  always 
vanishes.  Equations  (3*12)  and  (3-13)  reduce  to 


B22  U’xxxx"  £B26 
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Another  similar  special  case  when  all  the  pairs  are  running  either  parallel 
or  perpendicular  to  the  coordinate  axes  permits  simplification  of  the  equations . 
The  elastic  coefficients  in  this  case  are 
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9  in  this  case  is  either  0°  or  90°.  All  the  A.^  and  survive  except  the 
ones  with  subscripts  l£  and  26.  Consequently  ana  a^°  vanish.  The 
governing  equations  here  are 
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These  equations  describe  a  plate  which  is  orthotropic  in  both  in-plane  and  bending 
effects.  Furthermore,  if  there  are  the  same  number  of  pairs  in  both  directions, 
then  An  “  A22  and  B11  =  ®22;  80  Eq,  (4. 28)  becomes  the  same  as  Eq  (4.26). 

Several  authors  have  dealt  with  equations  of  the  same  form  as  Eqs  (4.24) 
and  (4.25)  in  connection  with  particular  problems.  Green  (ll)  exhibited  the 
method  of  solution  to  these  equations  in  connection  with  an  aeolotropic  single- 
layer  plate  subjected  to  in-plane  forces.  Lekhnitskii  (12),  (13)  also  solved 
the  same  pair  of  equations .  Both  these  authors  used  a  technique  involving  complex 
variable  theory.  Pell  (l4)  solved  a  problem  of  thermal  stresses  in  a  thin  aniso¬ 
tropic  plate.  Luxenberg  (15)  solved  a  problem  of  the  torsion  of  an  anisotropic 
plate.  Both  of  these  problems  were  governed  by  an  equation  of  the  same  form  as 
Eq  (4.25).  Other  authors,  for  example  Green  and  Taylor  (16)  (17),  Green  (18,  19), 
Fridman  (20),  Okubu  (21),  Morris  (22),  have  dealt  with  equations  of  the  form  of 
Eqs  (4.28)  and  (4.29).  They  essentially  used  the  same  method  as  that  used  by 
Green  and  Lekhnitskii.  Girkmann  (23)  discusses  the  solution  to  Eq  (4.29)  using 
a  double  Fourier  Series  in  connection  with  an  orthotropic  single-layer  plate; 
his  method  will  be  discussed  in  Section  7  of  this  report. 

When  the  equations  are  of  the  form  of  Eqs  (4.24)  and  (4.25),  the  distribution 
of  stress  near  a  hole  in  an  infinite  laminated  anisotropic  plate  can  be  calculated 
using  the  technique  presented  by  Savin  (24)  or,  for  the  case  of  a  circular  hole 
in  a  laminated  plate  which  behaves  orthotropically,  those  presented  by  Green  and 
Taylor(l7).  Savin  extended  the  Muskhelishvili  method  for  the  solution  of  plane 
elasticity  problems  to  solve  the  anisotropic  plane  problem  for  an  infinite  plate. 
Since  in  the  absense  of  transverse  loads,  the  differential  equations  for  the  plate 
and  for  plane  elasticity  are  of  the  same  foirn,  the  method  is  applicable  to  both 
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cases.  The  determination  of  general  expressions  for  stresses  is  not  practical 
since  the  stresses  are  determined  from  only  the  real  part  of  the  complex  poten¬ 
tials  which  comprise  the  solution. 


SPIRAL  ORIENTATION 

This  classification  refers  to  a  layered  plate  made  of  laminas  such  that  the 
elastic  axes  are  distributed  in  the  form  of  a  spiral.  Only  the  case  denoted  as 
equiangular  spiral  orientation  will  be  discussed  since  all  other  cases  are  too 
complicated  to  be  investigated  on  a  general  level.  The  equiangular  spiral 
orientation  describes  a  plate  formed  by  alternating  succeeding  laminas  from  the 
middle  surface  at  equal  but  opposite  angles  (one  positive  and  one  negative) 
frcm  a  coordinate  axis.  Each  succeeding  pair  must  have  a  larger  angle  than  the 
previous  pair.  For  this  case  all  the  Ay  and  Dy  survive  except  those  with 
subscripts  1 6  and  26.  All  the  Dy  are  annihilated  except  for  the  ones  with 
subscripts  l6  and  26.  The  determinative  equations  become 
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These  equations  are  partly  coupled  since  both  dependent  variables  appear  in  the 
equations . 

The  solutions  to  Eqs  (4.30)  and  (4.31)  nay  be  found  by  raising  the  order 
of  the  differential  equations  up  to  the  eighth  order  as  was  discussed  in  the 
beginning  of  this  section. 


SPECIALIZATION  OF  THE  CYLINDRICAL  SHELL  EQUATIONS  FOR  PARTICULAR  ORIENTATIONS 

The  two  determinative  Eqs  (2.l£)  and  (2.17)  can  be  specialized  for  certain 
orientations  of  the  laminas  composing  the  layered  shell.  The  following  discussion 
will  be  restricted  to  cylindrical  shells  composed  of  laminas  of  the  same  thickness. 
Suppose  that  the  laminas  are  situated  such  that  one  lamina  of  a  pair  is  at  the 
same  distance  above  as  the  other  is  below  the  middle  surface  of  the  shell  (now 
taken  as  the  reference  surface).  Both  laminas  have  the  same  orientation  of  their 
principal  elastic  axes.  In  this  instance,  the  Dy  always  vanish.  Accordingly, 

Eqs  (2.16)  and  (2.17)  reduce  to 
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(4.33) 
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In  the  above  we  have  set 

i  5  =  a 

a  elf  “  oy 


These  equations  indicate  that  the  laminated  shell  is  equivalent  to  a  homogeneous 
anisotropic  shell.  Further  simplification  is  possible  by  having  two  pairs  of 
Inminas  situated  so  that  there  exists  symmetry  about  the  coordinate  axes  x,  y. 

In  this  case  (called  pairwise  equiangular  orientation)  A^jg  ■  “  0  ,  thus 

Big  =  Bgg  =  0  in  Eq  (4.32).  A  further  specialization  results  if  all  pairs  are 
situated  at  plus  or  minus  45  degrees  from  the  coordinate  axes.  In  this  instance 
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Another  special  case  arises  when  all  pairs  are  situated  either  parallel  or 
perpendicular  to  the  coordinate  axes.  The  governing  equations  here  take  the  form 
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These  equations  describe  a  laminated  cylindrical  shell  which  is  orthotropic  in 
both  in-plane  and  bending  effects. 

Equations  (4.32)  and  (4.33),  or  the  specializations  discussed  subsequently, 
are  of  the  general  form 
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where  L-,  and  Lq  are  fourth  order  linear  differential  operators  in  x,  y.  Special 
forms  of  the  operators  occur  for  particular  orientations  of  the  1  ami  nan  as 
described.  A  method  of  solving  equations  of  this  type  has  been  suggested  by 


47 
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In  the  above  we  have  set 
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These  equations  indicate  that  the  laminated  shell  is  equivalent  to  a  homogeneous 
anisotropic  shell.  Further  simplification  is  possible  by  having  two  pairs  of 
laminas  situated  so  that  there  exists  symmetry  about  the  coordinate  axes  x,  y. 

In  this  case  (called  pairwise  equiangular  orientation)  A^g  =  m  0  ,  thus 
®l6  “  ®26  =  0  in  Eq  ( ^4-.  32) .  A  further  specialization  results  if  all  pairs  are 
situated  at  plus  or  minus  h-5  degrees  from  the  coordinate  axes.  In  this  instance 
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Another  special  case  arises  when  all  pairs  are  situated  either  parallel  or 
perpendicular  to  the  coordinate  axes.  The  governing  equations  here  take  the  form 
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These  equations  describe  a  laminated  cylindrical  shell  which  is  orthotropic  in 
both  in-plane  and  bending  effects. 

Equations  (4.32)  and  (4.33),  or  the  specializations  discussed  subsequently, 
are  of  the  general  form 
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where  Ln  and  Lg  are  fourth  order  linear  differential  operators  in  x,  y.  Special 
forms  of  the  operators  occur  for  particular  orientations  of  the  laminas  as 
described.  A  method  of  solving  equations  of  this  type  has  been  suggested  by 
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Vlasov  (7)-  Assume  the  existence  cf  a  function 


y)  defined  such  that 
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then  Eq  (1.33)  is  substituted  into  Eqs  (1.  30)  and  (1.37),  the  first  equation  is 
identically  satisfied,  and  the  second  equation  tahes  the  form 
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Stress-resultants,  stress-couples  or  reference  surface  strains  can  be  expressed 
in  terns  of  ^  in  a  straight-forward  manner. 

For  the  case  of  a  laminated  orthctrcpic  cylindrical  shell,  Eqs  (I.3I)  and 
(1. 35 )  are  replaced  by  Eq  (1.39)  vliere  L,,  L-,  take  the  fora 


L1  *  B22  ^  ^’xxxx  +^B66  +  2B12^  ^  ^xxyy  +  B11  ^  ^'yyyy 


Lo  =  Dpp  (  ),yaggx  +  2(D10  +  2D,^)  (  ),_  +  D„  (  ), 


12  06 J 


xxyy 


yyyy 


Finally,  for  an  isotropic  shell, 
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where  E,  3^  are  modulus  of  elasticity  and  Poisson's  ratio  and  h  the  plate 
thickness . 
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PART  II:  APPLICATIONS  OF  PLATE  THEORY 


INTRODUCTION 

Several  applications  of  the  general  theory  for  laminated  anisotropic  plates 
are  discussed  in  this  part.  First,  the  equations  of  the  plate  theory  are  modified 
for  the  cylindrical  bending  and  extension  of  long  rectangular  plates  so  that  solu¬ 
tions  may  be  obtained  by  integration  or  by  transposition  of  known  solutions  from 
beam  theory.  One  example  is  given  to  illustrate  a  step  by  step  method  of  obtain¬ 
ing  the  elastic  constants  for  a  plate.  Next,  solutions  are  obtained  by  determin¬ 
ing  the  edge  loads  which  will  cause  the  deflected  surface  to  become  a  prescribed 
quadratic  surface.  The  following  examples  sire  given:  uniform  tension,  pure 
shear,  uniform  bending  about  one  axis,  and  pure  twist.  As  discussed  under  methods 
of  solution,  for  certain  orientations  the  general  equations  reduce  to  those  of 
orthotropic  plate  theory.  The  solution  for  the  uniformly- loaded  s iraply-supported 
rectangular  plate  is  presented.  As  illustrations  of  the  solution  of  general 
equations  by  the  perturbation  method,  three  examples  are  given:  the  clamped 
circular  plate  (for  which  an  exact  solution  can  also  be  obtained),  the  simply- 
supported  rectangular  plate,  and  the  rectangular  plate  simply  supported  on  two 
opposite  edges  and  clamped  on  the  other  two  edges.  Finally,  the  problem  of 
obtaining  the  optimum  arrangement  and  orientation  of  the  layers  so  as  to  minimize 
the  deflection  in  rectangular  plates  is  discussed. 


5-  WEAK  CYLIinjRICAL  BEEDUiG  AED  EXTENSION  OF  A  LONG 
RECTANGULAR  LAMINATED  FLATS  OF  ORTHOTROPIC  MATERIALS 


A  class  of  problems  vhich  deserves  attention  is  that  of  extensional  and 
flexural  behavior  of  a  long  rectangular  laminated  plate  composed  of  orthotropic 
material,  When  such  a  plate  is  loaded  by  a  system  of  forces  which  does  not  vary 
in  the  y-direction,  the  deflected  surface  at  same  distance  from  the  ends  is 
cylindrical  (i.e.,  it  is  also  independent  of  y).  Therefore,  the  displacement 
v  is  zero  and  the  displacements  u  and  w  are  independent  of  y.  The  plate  may  be 
investigated  by  studying  a  typical  elemental  strip  bounded  by  tvo  planes  y  = 
const.  (Figure  5.1).  In  essence  the  problem  becomes  one-dimensional,  and  the 
deflection  curve  resembles  that  of  a  deflected  beam. 


As  a  consequence  of  the  type  of  displacements,  6  y  =  0  and  =  0.  The 

stress-strain  relations  given  by  Eq  ( 1. 6)  simplify  considerably  for  the  elemental 
plate  strip. 
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For  this  plate  strip,  the  differential  equations  for  the  transverse  deflec 
tion  w  and.  the  in-plane  displacement  Uq  can  be  uncoupled.  Two  general  cases  of 
loading  will  be  discussed:  weak  cylindrical  bending  and  uniaxial  extension. 


WEAK  CYLINDRICAL  BENDING  OF  A  PLATE  STRIP 

Consider  a  simply  supported  plate  strip  of  length  L.  The  origin  of  the 
coordinate  system  is  located  at  the  left  hand  side  of  the  plate  strip.  The  x-y 
plane  lies  in  the  reference  surface  of  the  plate;  the  x  and  z  axes  are  directed 
to  the  right  and  downward,  respectively,  as  shown  in  Figure  5.2.  This  plate 
strip  is  loaded  by  an  arbitrary  system  of  transverse  forces  and  moments  M-^  and 
Mg  at  the  ends.  The  in-plane  force  K*  is  taken  to  be  zero. 
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Under  the  action  of  the  external  loads,  the  plate  strip  deflects.  The 
deflection  curve  can  be  determined  in  the  following  manner.  Since  1CX  is  cere, 
by  setting  the  first  of  Eq  (5-2)  equal  to  cero,  the  following  relation  between 
Uq  and  w  is  obtained. 
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Substitution  of  the  above  equation  into  the  first  of  3q  (5-5)  gives 
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Equation  (5-5)  is  the  governing  differential  equation  for  the  deflection  curve 
w  and  Sq  (5-M  is  the  associated  differential  equation  for  the  in-plane  reference 
surface  displacement  uQ.  Equation  (5-5)  is  of  the  same  form  as  the  equation  for 


the  deflection  of  an  isotropic  beam. 

(°ll)2  _  A11  D11 


Only  the  flexural  rigidities  are  different. 


The  quantity 
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replaces  the  usual  flexural  rigidity 


cf  the 


isotropic  beam.  Therefore,  transposition  of  known  solutions  for  isotropic  beams 
with  the  corresponding  loading  conditions  is  possible.  Hie  solution  may  also  be 
obtained  by  double  integration  if  the  expression  for  Mj_(x)  is  known. 

The  solution  to  Eq  (5 • M  may  be  obtained  by  direct  integration 
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The  constant  k  is  to  be  evaluated  from  a  prescribed  boundary  condition.  The 
boundary  condition  for  u^,  can,  however,  be  arbitrarily  selected  since  it  repre¬ 
sents  only  a  point  from  which  the  displacement  is  measured.  Only  the  derivative 
UqjX  which  is  the  in-plane  strain  is  of  any  interest  in  this  problem.  Therefore, 
let  the  boundary  condition  be  taken  as 
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Thus  far  the  problem  has  proceeded  as  for  the  case  of  the  cylindrical  bend¬ 
ing  of  a  single-layer  isotropic  plate  strip  except  for  the  different  flexural 
rigidities.  Now,  however,  unlike  the  single-layer  isotropic  plate  strip,  which 
requires  only  a  transverse  bending  moment  My.  to  maintain  the  cylindrical  deflec¬ 
tion  surface  during  bending,  additional  force  and  moment  components  (K^,  Ny, 
IJxy)  must  be  present  to  preserve  the  cylindrical  surface  for  the  laminated  J 
anisotropic  plate  strip.  Unless  these  forces  and  moments  are  developed  at  the 
supports,  the  deflected  surface  will  not  remain  cylindrical  and  warpage  will 
occur  along  the  y-di recoion.  These  forces  and  moments  may  be  computed  from 
3qs  (p.2)  and  (5- 3)*  The  stresses  may  be  computed  from  2q  (3*l3a). 

UNIAXIAL  3XT3NSI0N 

Consider  the  same  simply  supported  plate  strip  with  the  coordinate  system 
as  shown  in  Figure  5 • 3-  At  the  ends  of  this  plate  strip  an  in-plane  extensional 
force  of  magnitude  N  is  applied.  The  bending  moment  is  taken  to  be  zero 
in  this  case. 


Figure  5-3 


3ecause  this  force  is  applied  at  the  reference  surface  (which  is  also  the 
surface  of  the  coordinate  plane  x-y),  the  plate  strip  will  not  only  elongate  but 
also  may  bend.  This  behavior  occurs  because  the  applied  force  may  not  coincide 
with  the  centroid  for  a  tensile  or  ccmpressive  force  (the  point  at  which  an  axial 
force  causes  no  bending)  on  thi3  section.  The  deflection  curve  and  the  extension 
can  be  determined  in  a  manner  analogous  to  that  used  in  the  previous  case.  The 
first  of  3q  (5-3)  is  set  equal  to  zero,  from  which 


o,x 


(5-12) 


Substituting  3q  (5.12)  into  the  first  of  3q  (5*2)  gives  the  differential  equation 
of  the  deflection  curve  w. 
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"1 1'  J 


D 

.All  °11  " 


Again  Eq  (5*13)  is  of  the  same  form  as  the  equation  for  the  deflection  of 

are  different.  By  replacing 

for  the  anisotropic  laminated 
plate,  knovn  solutions  from  beam  theory  for  corresponding  loadings  may  be  used. 


an  isotropic  beam.  Only  the  flexural  rigidities 

*\2 


El  of  the  beam  equation  with 


A11  D11  "  (Dli)‘ 
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The  solution  to  Eq  (5-12)  is  obtained  by  direct  integration: 


4  k 


(5-14) 


Let  the  boundary  condition  for  the  evaluation  of  constant  k  be  arbitrarily 
selected  sis 

uo(°)  -  -  \  (5.15) 

The  reason  for  such  a  specification  will  be  apparent  shortly.  Therefore  k 
is  given  by 


k  *  k1  (5.16) 

To  maintain  a  flat  surface  during  extension,  a  bending  moment  must  be  applied 
to  cancel  the  deflection  caused  by  the  in-plane  extenslonal  force  acting  at  the 
reference  surface.  Application  of  this  bending  moment  may  be  thought  of  as  the 
transfer  of  the  applied  force  from  the  reference  surface  to  the  centroid  for 
tension-compression  of  this  section  as  shown  in  Figure  5.4. 


rfcLt-v 


O'. 


U 


centroid  for 
tensile  or 

oppressive  for  eg,  ~U£ 


Figure  5.4 


To  remove  the  deflection  curve,  a  moment 

at  the  ends.  Superposing  both  solutions  of  w 
plate. 


of  magnitude 


N  D 


’ll 


All 


is  applied 


then  gives  the  condition  of  a  flat 


v(x)  =  v  ,  (x)  +  v,  ,(:c)  =  0 
'  ext  '  bend 


and 


(x)  = 


bend 


D11  N 


(Aii  [(Dii)2  '  An  Dn] 

The  boundary  condition  for  uo-oen(i  is  taken  as 


x  +  k 


(5.17) 


(5.13) 


u  (o)  =  k 
°bend 


(5.19) 


so  that  for  the  complete  problem  the  displacement  %  at  x  =  o  is  zero.  Therefore 


"^bend 


(5.20) 


and 


o. 


(x) 


bend 


11 


N  (D-^) 
'(Dil)2  - 


Dn  Au 


x  +  k 


(5-21) 


The  superposition  of  u0(x)  for  extension  and  u  (x)  for  bending  gives 


u  (x)  =  ——  x 
o 


(5-22) 


As  was  mentioned  previously  the  transfer  of  the  in-plane  extensional  force  to 
the  centroid  for  tension-compression  of  this  particular  section  gives  the  same 
result.  The  eccentricity  of  this  centroid  from  the  reference  surface  is  given 
by  the  formula 


ecc 


N 


x 


M 

N 


11 

11 


(5-23) 


When  the  load  is  applied  at  the  centroid,  no  bending  will  occur. 


As  for  the  previous  case  of  cylindrical  bending,  additional  restraining  forces 
and  moments  must  be  developed  at  the  supports  to  prevent  varpage  along  the  y- 
direction.  IKiese  quantities  are  again  computed  from  Eqs  (5*2)  and  (5-3)  Stresses 
are  given  by  Eq  ( 3  •  l8a )  • 

'Cable  I  is  a  summary  of  the  cylindrical  bending  and  uniaxial  extension 
problems . 
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TABLE  1:  Summary  of  Cylindrical  Bending  and  Uniaxial  Extension  Problems 
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Construction  -  Manual  of  the  A.I.S.C.,  Fifth  Ed.,  pp.  366  ff 


SUPERPOSITION  OF  SOLUTIONS 


Many  problems  with  different  loadings  and  different  support  conditions  may 
be  found  by  superposing  the  previous  cases.  As  an  example,  consider  the  following 
plate  strip  shown  in  Figure  5-5- 


z 


Figure  5-5 


This  problem  may  be  solved  by  superposing  the  cases  of  a  transversely  loaded 
plate  strip  and  bending  by  unequal  moments.  The  bending  moment  for  the  bending 
problem  is 


M  -  M 

Mx(x)  =  Mx  +  (  ■2--l-  ■  -  )  x 

The  bending  moment  for  the  triangularly  loaded  plate  strip  is 

«*<*>  - 

w(x)  for  both  cases  can  be  obtained  by  double  integration.  Then  these  two  solu¬ 
tions  must  be  superposed  in  such  a  way  that 

v,x(o)  =0  and  w,x(  L)  =  0 

These  two  conditions  determine  the  magnitudes  of  M1  and  M^. 

■  -  <«2  *  »i>  5  -  ir  ■ 0 
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V'X(L) 


M  L  ML 

+  -i-  +  — 


0 


Solving  these  two  equations  gives 


Ml  = 


kL' 

30 


M- 


Note  that  these  end  moments  are  the  same  as  that  for  the  isotropic  beam  of 
the  same  loading  and  boundary  condition.  Thus  the  solution  could  also  have  been 
obtained  by  transposition  of  the  isotropic  beam  solution. 


EXAMPLE  OF  A  PLATE  OF  12  ORTHOTROPIC  LAMINAS  WITH  SPIRAL  ORIENTATION 

Consider  a  12  layer  laminated  plate  of  spiral  orientation.  Each  lamina  has 
thickness  h  and  is  oriented  such  that  the  principal  elastic  axes  are  at  either 
+  45°  or  -  45°  with  the  coordinate  axes.  Let  the  reference  surface  be  between 
the  6th  and  7th  layers  as  shown  in  Figure  5-6. 


Layer  1  +45°  ' 

- 1 

2  -45° 

5h 

3  +45° 

4h 

4  -45° 

^  hl 

3h 

5  +45° 

- - 

- ■  T'vh 

k  6  -45° 

- 3 

n-r 

7  +^5°  [ 

— ih 

!  3  -45°  s 

;  2h 

9  +45°  ! 

3h 

10  -45°  « 

, 

4h 

11  +45°  ' 

h12 

!» 

12  -45°  ! 

6h 

t 

z 


Figure  5-6 


In  this  example  let  the  four  independent  parameters  be  E  ■  E,  k  -  1. 2, 

0.3,  and  X  =  0.42.  The  principal  elastic  constants  are  Cu  =  1^12 IE, 

C22  =  1. 345E,  C12  ■  0.4o4e  and  C66  =  0.471E.  The  elastic  constants  C±j  for  the 
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rotated  axes  are  computed  using  Eq  (l. 12). 


For  a 

lamina  at  +  45° 

For  a 

lamina  at  -  45° 

1.29  0.35  0.06 

1.29  0.35  -O.06 

*13  =E 

0.35  1.29  0.06 

*13  “E 

0.35  1.29  -0.06 

0.06  0.06  0.4l 

-0.06  -0.06  0.4l 

4 

*  /  V 

The  coefficients  A^.,  Dij,  and  are  computed  from  Eq  (1-30)  as  shown  in  Table  2 
and  may  be  summarized  as  rollows: 


M 


Eh 


15-48 

4.20 

0 

4.20 

15.48 

0 

0 

0 

4.92 

0 

0 

-0.36 

0 

0 

-0.36 

-0.36 

-0.36 

0 

'±3 


Eh~ 


185.76  50.44  0 

50.44  1B5.76  0 

0  0  59-08 


With  Aji,  Dji  4,  and  Dj<  known,  Bj4,  bjj,  T.,  ( 1  ^ ,  and  can  be  computed  from 

Eqs  (2.11),  (2.13),  (2.20),  and  (2.21).  J 


UiJ 


b 


ij 


=  h 


0.0697 

-0.0189 

0 

-0.0189 

0.0697 

0 

0 

0 

0.203 

0 

0 

-0.073 

0 

0 

-0.073 

-0.018 

-0.018 

0 

6l 


O.O83 

0.000+ 

±0.012 

m  U) 

ij 

1 

=  h 

0.000+ 

0.083 

+0.012 

±0.003 

±0.003 

0.083 

+0.001 

+0.001 

-0.119 

ij 

=  Eh 

+0.001 

+0.001 

-0. 119 

-0.007 

-0.007 

+0.009 

In  the  case  where  there  are  two  signs,  the  upper  sign  is  for  an  odd  numbered 
layer  and  the  lower  sign  is  for  an  even  numbered  layer. 

Three  cases  of  the  rectangular  laminated  plate  strip  for  this  spiral 
orientation  have  been  computed: 

(1)  Cylindrical  bending  by  a  uniform  moment. 

(2)  Uniaxial  extensions :  Since  the  eccentricity,  eccw  ,  is  zero 

x 

for  this  particular  location  of  the  reference  surface,  no  shifting  of 
the  applied  force  is  necessary. 

(3)  Uniformly  loaded  plate  strip. 

All  the  results  are  obtained  by  substituting  these  numerical  values  into  Table  1. 
The  final  results  are  summarized  in  Table  3-  Plot6  of  the  stress  distributions 
are  shown  in  Figure  5-7  and  Figure  5-8* 
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liable  3 :  Summary  of  the  Rectangular  Plate  Strip  Problem 
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6.  DEFLECTION  SURFACES  FOR  PRESCRIBED  EDGE  LOADS 


In  the  absence  of  transverse  loads  (i.e.,  q  =  0),  the  compatibility  and 
equilibrium  equations,  Eqs  (3.12)  and  (3*13)j  are  both  satisfied  identically  if 
U  and  w  are  assumed  as : 


2  2 

U  =  cqx  +  c^y  +  c2xy 


2  2 

w  =  gQX  +  g^y  +  ggXy 


(6.1) 


where  the  c^  are  coefficients  which  describe  the  form  of  the  stress  function  and 
the  gi  are  coefficients  which  describe  the  shape  of  the  deflection  surface. 

The  introduction  of  Eq  (6.1)  into  Eqs  (3-ll)  and  (3.10a)  leads  to  the  follow¬ 
ing  expressions  for  the  in-plane  forces  and  bending  moments: 


N  =  U,  =  2c, 
x  ’yy  l 


N  =  U,  =  2c 
y  xx  o 


(6.2) 


N  =-U,  =  -  c,, 

xy  ’xy  2 


x  —  1  “11  ““o  “12  "  C2  bl6  +  2go  ‘Hi  +  2gl  dI2  +  ^2  dl6 


M  =  2c,  b„  +  2c  b 
1  11  o 

M  =  2c.  b —  +  2c  b —  -  c_  b 


y  1  “21  ~“o  “22  “2  b26  +  ^o  d12  +  2gl  *22  +  2g2  d26 

Mxy  °  2cl  b6l  +  2co  b62  -  C2  b66  +  ^o  dl6  +  2gl  *26  +  2g2  d66 


(6.3) 


Substituting  the  values  of  the  ci  from  Eq  (6.2)  into  Eq  (6.3)  and  writing  the 
resulting  equation  in  matrix  fozra,  one  obtains: 


M  =  bN  +  2dg 


(6.4) 


where 


(6.5) 
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From  Eq  (6.k),  the  deflection  surface  coefficients  may  be  expressed  in  terns 
of  the  bending  moments  and  in-plane  forces. 


H 


M  -  bN 


(6.6) 


where  the  flexibility  matrix,  f,  is  related  to  the  stiffness  matrix,  d,  by: 


f  =  d 


(6.7) 


The  shape  of  the  deflection  surface  for  several  cases  of  edge  loadings  will 
now  be  considered. 

UNIFORM  TENSION 

Let  all  force  components  be  zero  except  3^,  then  the  g^  become: 


So  -  -  5  [  fli  bll  +  b21  +  fl£  b6i] 

®1  ■  ■  I  "*  [  f12  bll  +  f22  b21  *  f2«  b6l] 

e2  *  "  2  Hx  £  fJjS  bll  +  f26  b21  +  f66  b6lj 


(6.8) 


and  the  deflection  surface  is  given  by: 


'  -  -  5  »,  +  f12»21  *  fl6b6l>  +  <f12bU  +  f22b21  +  r26b6l>  * 

*  (fl£bll  +  f26b21  +  f66b6l>  xy]  <6'9) 


It  must  be  recalled  that  all  forces  are  applied  at  the  reference  surface  and 
that  it  may  be  possible  to  change  the  fora  of  the  deflection  surface  by  changing 
the  location  of  the  reference  surface;  that  is,  one  of  the  coefficients  of  x  , 
y^,  or  xy  may  be  made  to  vanish  by  a  particular  choice  of  the  reference  surface. 

PURE  SHEAR 

Assume  that  the  only  non- zero  force  component  is  N^.  The  displacement 
surface  coefficients  are  given  by: 


go  ”  "  2  Nxy  £fll  bl6  +  f12  b26  +  fl6  to66  J 
S1  =  ■  2  Nxy  £f12  bl6  +  f22  b26  +  f26  ^66^ 


(6.10) 
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(6.io) 


g. 


s2  2 
and  the  displacement  surface  is 


fl6  bl6  +  f26  b26 


+  f66  b66 j 


2 
p 

(6.11) 


»  ■  -  i  v  [<rilbl6  +  VaS  +  fl6bM>  +  <fV,16  +  f22b26  +  f26  b66>  * 

+(fl6blS  +  r26  b26  +  f66  b66>  **] 

As  in  the  previous  case,  one  of  the  coefficients  of  x^,  y2,  or  xy  may  be  made 
to  vanish  by  an  appropriate  choice  of  reference  surface. 

UNIFORM  BENDING  ABOUT  ONE  AXIS 

Let  all  of  the  force  components  be  zero  except  which  is  assumed  to  be 
constant.  Then,  the  values  of  the  g^  are  given  by: 

«o  ■  I  fll  M* 


«1  ■  5  f12  Mx 


(6.12) 


*2  =  2  f 16  Mx 


and  the  deflection  surface  becomes : 

w  =  i  MX  [  fil  x2  +  fl2  y2  +  fi6  ]  (6.13) 

In  this  case,  the  only  tern  which  can  possibly  be  made  to  vanish  by  changing 
the  location  of  the  reference  surface  is  fj£.  The  other  terms  must  exist  for  all 
possible  choices  of  the  reference  surface  and  for  all  possible  orientations  of  the 
laminas. 


PURE  TWIST 


Assume  that  the 
the  form: 


only  force  component  which  exists  is 


then  the  take  on 


6o  "  2  fl6  Mxy 
S1  =  2  f26  Mxy 


(6.14) 
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(6.14) 


e2  2  f66  Mxy 

The  deflection  surface  is  given  by 

w  =  i  Mxy  [  flS  x2  +  f26  y2  +  f66 

The  only  terns  which  might  be  made  to  vanish  in  this  case  are  fj£  and/or  f^. 
Both  the  location  of  the  reference  surface  and  the  orientation  of  the  laiUnaB  will 
affect  the  existence  of  these  terns. 
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7.  BENDING  OF  ORTHOTROPIC  LAMINATED  PLATES 

As  remarked  in  Section  4,  the  general  determinative  Eqs  (3-12)  and  (3.13) 
governing  the  behavior  of  laminated  plates  can  be  specialized  in  certain  instances 
of  lamina  orientation.  In  the  particular  case  of  a  plate  with  pairwise  orientation 
whose  layers  are  90  degrees  relative  to  one  another,  the  in-plane  and  bending 
effects  are  uncoupled.  In  this  case  a  solution  of  Eq  (3*13)>  now  in  the  special 
form  of  Eq  (4.29),  may  be  obtained  by  the  Navier  method  (see  for  example.  Reference 
(23)).  Recalling  the  form  of  this  equation: 

D,_  w,  +  (2D. _  +  4D£/;)  w,  +  D_0  w,  =  q(x,y)  (7-l) 

11  ’xxxx  12  66'  'xxyy  22  yyyy  v  ' 

if  the  load  q(x,  y)  is  ejtpanded  in  a  double  Fourier  series 


q(x,y) 


y 


f  a  sin^ 


(7.2) 


then  the  solution  form  of  w(x,y)  for  the  simply  supported  plate  may  also  be  taken 
as  a  double  Fourier  series 


v(x,y)  .  N  V  w  smiBi  si„S££ 
'  nm  a  b 

n=l  m=l 


(7.3) 


where  a  and  b  are  the  lengths  of  the  sides  of  the  plate.  The  unknown  coefficients 
wmn  may  be  found  by  substituting  the  two  series  back  into  Eq  (7-1)  and  equating 
the  like  terms.  The  equation  is  satisfied  for  every  value  of  n  and  m  if 


nm 


D 

* 


*7 


(2D 


12 


“V 


2  2 
n  m 

2  2 
a  b 


IT 

n  m 
D22  ^ 
D 


(7.4) 


The  complete  solution  of  w(x, y)  is  then, 


.<*,>. T 


n=l  m=l 


q  sin  2^  sin 

Tim _ a _ b _ 

D11  T  +  2<D12  +  2D66>  %  + 

a  abb 


(7.5) 


for  the  special  case  of  a  uniform  load  q  =  q^  q^  is  given  by 
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(7-6) 


m  =  1)  3)  5,  -  -  - 


n  =  1,  3,  5, - 

The  moments  and  stresses  are  calculated  from  Eqs  (3.1Da)  and  (3.18a).  In¬ 
plane  forces  which  are  found  from  the  solution  to  Eq  (4.28)  may  he  superimposed 
as  a  separate  problem  since  the  differential  equations  and  boundary  conditions 
are  uncoupled  for  this  particular  orientation  of  the  laminas. 

EXAMPLE  OF  A  SQUARE  PLATE  OF  13  LAYERS  WITH  PAIRWISE 
90  DEGREE  EQUIANGULAR  ORIENTATION 

Consider  a  13  layer  laminated  plate  of  such  an  orientation  that  its  behavior 
can  be  described  by  Eqs  (4.28)  and  (4.29).  Each  lamina  has  thickness  h  and  is 
orientated  such  that  the  principal  elastic  axes  coincide  with  the  coordinate  axes. 
The  length  of  the  sides  is  a.  Let  the  reference  surface  be  the  middle  surface 
of  the  seventh  layer. 

-6.5h 


-►  x 


6.5h 

■  I 

I 

T 

z 

Figure  7.1 
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Snn  =  ~~2  ~ 

n  mn 


Let  the  four  independent  parameters  be  the  same  as  that  of  the  example  in 
Section  5,  E  “  E,  k  ■  1. 2,  V  “0.3;  and.  \  ■  0.42.  The  are  given  by  Eq  (l.  14). 
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for  an  even  numbered  layer 


for  an  odd  numbered  layer 


► 

1. 345 

O.4o4 

0 

1.121 

0.404 

0 

0.4o4 

1. 121 

0 

ciJ =  E 

o.4o4 

1.345 

0 

0 

» 

0 

0.471 

0 

0 

0.471 

The  coefficients  Ajj  and  are  given  by  Eq  (l.  30).  is  zero  for  this 

orientation.  A  convenient  form  such  as  Table  2  of  Section  5  may  be  used  for  the 
computation  of  these  coefficients.  The  results  are 


A.  .  =  Eh 
ij 


15.917  5.252  0 

5.252  l£.l4l  0 

0  0  6.123 


221. 028 
73.966 
0 


73.966 

230.455 

0 


0 

0 

36.232 


Knowing  A^j,  the  Bj^  may  be  found. 


0.0704 


0.0229 


0 


0.0229  0 

0.0694  0 

0  0. 1633 


The  governing  differential  Eqs  (4.28)  and  (4.29)  became 

~  T  0.0694  U,  +  0.1175  U,  +  0.0704  U,  1=0 

Eh  I  xxxx  xxyy  yyyy  J 

Eh3  221.028  w,  +  492.860  w,  +  230.455  w,  1  =  q(x,y) 

L  xxxx  '  xxyy  yyyy 

Since  these  equations  are  uncoupled,  their  solutions  may  be  taken  independent¬ 
ly.  Only  the  bending  problem  will  be  considered  here.  If  the  plate  has  a  uniform 
load  q  =  qQ,  then  the  solution  for  w(x,y),  which  is  given  by  Eq  (7*5)  and  Eq  (7-6) 
is 
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v(x,  y) 


16  q.  V  V  sin  2£sin  5SJT 

o  /  /  _ a _ a _ 

jt6a4Eh3  n=l  m<L  221.028  n5m  +  492.860  m3n3  +  230.455  mP 

(n  -  1,  3,  5 - ;  o  «  1,  3,  5 - ) 

The  moments  sure  given  by  Eq  (3. 10a) 

Mx  -  -  (221.028  v,^  +  73.966  v,  )  Eh3 

My  =  -  (73.966  +  230.455  V>yy)  ^ 

M  -  -  ( 172. 464  w,  \  Eh3 
xy  V  'xy  / 

The  stresses  are  given  by  Eq  (3-lSa).  Since  the  in-plane  forces  and  t^^ 
are  zero,  these  equations  reduce  to 


a  W  =  -  z 

X 

'  p  (i)  „  +  7;  (i)  „  ] 

_  C11  V’xx  +  C12  w’yyJ 

a  (i)  -  -  z 

y 

C  w  +  'c  w  1 

_  12  ’xx  °22  'yyl 

C\J 

l 

n 

•H 

C  (1)  v 

w,.  . 

&  & 

The  deflections  and  the  moments  at  the  center  of  the  plate,  x  =  y  =  — , 
are  summarized  below 


terns  in  series 

(nm) 

w 

M 

X 

M 

y 

11 

4 

q  a 

0.0000176  * 

Eh-5 

0.0513  qQ  a2 

0.0530  qo  a2 

11,  13,  31,  33, 

4 

q  a 

0.0000172  -2— 

0.0466  qQ  a2 

0.0481  qQ  a2 

51,  15 

EhJ 

Additional  terms  are  needed  to  obtain  more  exact  values  for  the  moments 
because  their  series  representations  converge  slower  than  the  series  representa¬ 
tion  of  the  deflection. 
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for  odd  numbered  layer 


)  (0.0000766) 


for  even  numbered  layer 


8.  APPLICATIONS  OF  THE  PERTURBATION  METHOD  OF  SOLUTION 

It  should  be  recalled  that  the  perturbation  method  of  solution,  as  described 
in  Section  4,  is  applicable  when  the  squares,  products,  and  higher  powers  of  the 
perturbation  parameters  are  negligibly  small. 

The  perturbation  method  of  solution  may  be  applied  to  problems  in  which  the 
boundaries  are  not  rectangular  with  little  more  difficulty  than  arises  in  applying 
it  to  problems  in  which  the  boundaries  are  rectangular.  This  situation  arises 
because  the  differential  equations  for  the  first  approximation  and  for  each 
correction  tem  are  of  the  same  form  as  the  equations  for  isotropic  plates.  Hence, 
any  problem  which  has  been  solved  for  isotropic  plates  can  be  extended  to  laminated 
plates  by  the  perturbation  method.  As  an  example  of  a  problem  for  which  the 
boundaries  are  not  rectangular,  the  solution  for  a  uniformly- loaded  laminated 
circular  plate  with  clamped  edge  is  indicated  below.  An  exact  solution  can  be 
obtained  for  this  case,  so  that  a  comparison  can  be  made  between  the  exact  solu¬ 
tion  and  the  perturbation  solution.  In  addition  to  this  solution,  the  perturbation 
solutions  for  rectangular  plates  simply  supported  on  all  four  edges  or  simply 
supported  on  two  opposite  edges  and  clamped  on  the  other  two  edges  are  also 
indicated. 


CLAMPED  CIRCULAR  PLATE 


From  Sq  (4.15),  the  differential  equations  for  the  first  approximation  are 
the  same  as  those  for  isotropic  plates 


v  o  D 

4 

TT  U  =  0 

v  o 


(8.1) 


The  boundary  conditions  for  a  clamped  edge  will  be  taken  to  be: 


wl 


'  r 


r=a 


r=a 


=  0 


-  U,  +  —  U, 
r  r  ’ee 


r=a 


N 


=  0 


r=a 


For  a  uniform  load,  the  solutions  of  these  differential  equations  in  polar 
coordinates  are:  (see,  for  example.  Reference  (25;  ) 


w  =  ^o  /  2  2s2 

°  555  (a  ‘  r  > 


u  =  o 

O 


(8.2) 
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where  qQ  is  the  magnitude  of  the  load  and  a  is  the  radius  of  the  plate. 

When  Eqs  (4.l6),  (4.17),  (4.18),  and  (4.19)  are  transformed  into  polar 
coordinates,  the  following  equations  are  obtained: 


4 

V 

uill 

=  *  Ui22 

=  -«** 

55 

v4 

u*66 

=  *0 

55 

V 

U412 

4  „ 

=  *  Ui21 

-  -  3qo 

55” 

v4 

Wnll 

4 

=  V  wnl!  1 

=  3<1o 

55” 

4 

V 

w  * 
nl2 

=  «o 

55 

V4 

wn66 

=  v4  wn6g  = 

(1  -  vD) 

v4 

Wn22 

“  V  wn22  = 

3*o 

55” 

(8.3) 


(8.4) 


The  remaining  equations  are  all  of  the  form  v  (U,v)  =  0,  with  their  solutions 
identically  equal  to  zero.  The  solutions  of  the  above  equations  are: 


UAll 

ui66 


U 


122 


%  ,  2  2x2 

512D  (a  -  r  ) 


%  ,  2  2x2 

(a  -  r  ) 


(8.5) 


uii2  " 

Ui21 

»  .  £5°  (a: 

512D  '■a 

3<10  (a2 

=  512D  ^ 

wnll  = 

W  -  T 

nil 

TJ  ^  —5 

*0 

CM 

y-’s 

1 

Wnl2 

512D 

wn 66  " 

wn<£ 

-  <x  -  VB>  ' 

2  2x2 


2x2 


(8.6) 


2x2 
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„  J*0  /  2  2x2 

Wn22  =  Wn22  “  512D  "  r  ) 


(8.6) 


Therefore,  the  complete  second  approximation  is : 

rrmm  (&2  -  2)2  [  1  +  (nn  +  nu  +  n22  +  n22>  \ 

1  1  -*D 

+  nil  5  +  (n66  +  n6&  ( — ET-^ 


(8.7) 


qo  ,  2  2.2 

(a  -  r  ) 


£  “^66  ”  +  ^22}  2  "  ^12  +  ^1^  2  ]  (8.8) 


An  exact  solution  for  the  uniformly- loaded  clamped  circular  plate  may  be 
obtained  by  recognizing  that  the  following  forms  of  solution  satisfy  boundary 
conditions : 


„  /  2  2x2 

v  =  C1  (a  -  r  ) 


TT  „  /  2  2x2 

U  =  (a  -  r  ) 


(8.9) 


Two  simultaneous  equations  for  the  unknown  constants  C-,  and  Cp  are  obtained  by 
substituting  the  above  forms  of  solution  into  Eqs  (3. 12)  and  (3*  13).  The  following 
solutions  are  obtained  when  the  constants  C-^  and  C2  have  been  determined: 


U  -  *2  (ac  -  r  ) 


(a2  -  r2)2  [3B22  +  ^12  +  B66>  +  3Bllj 
0  ,_2  _  r2j2  pb12  *  ^bll  +  b22  2b66^  *  3b2lj 


(8.10) 


where. 


13*22  +  (2*12  +  *66>  +  3*11  ][  3dU  +  2<412 


3*11  +  2(*„  + 


2466>  +  3*22] 


3*12  *  (*11  +  *22  -  2b66)  +  3*2 


When  this  exact  solution  is  expressed  in  terms  of  the  perturbation  parameters 
and  the  division  is  expressed  in  series  form,  the  perturbation  solution  can  be 
seen  to  consist  of  the  first  terms  of  the  exact  solution. 
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The  solution  for  the  simply  supported  circular  plate  can  be  obtained  in  a 
similar  manner  provided  due  care  is  taken  in  properly  satisfying  the  boundary 
conditions . 

SIMPLY-SUPPORTED  RECTANGULAR  PLATE 

For  the  present,  the  load  will  be  assumed  to  be  sinusoidally  distributed; 
that  is 


-  %  sin  r  sin  ? 


(8.11) 


where  the  coordinate  system  is  as  shown: 


x 


b 


a 


y 


This  will  serve  to  illustrate  the  method,  and  will  indicate  that  the  solution 
for  any  load  distribution  may  be  obtained  by  expanding  the  load  in  a  double 
Fourier  series.  For  this  load  distribution,  the  solution  to  the  first  approxi¬ 
mations,  i.e.,  the  isotropic  case  is:  (see  Reference  (25),  for  example) 

a4 

w  ~  v  1  q-°' - 2~~2  sin  sin  &  (8.12) 

°  «*  D  <1  +  a2)2  *  h 

„o=  0 


where  a  *  a/b,  the  aspect  ratio. 

Inserting  these  relationships  into  Eqs  (4.l6),  (4.17),  (4.18),  and  (4.19), 
the  following  equations  are  obtained: 
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Vk  Vl  -  V*  Ofc* 


% 


ac 


D  (l  +  a  ) 


— Sin  E2S  Sin  £Z 

A  \  &  D 


\7  u. 


'466 


2q  a‘ 

— °— g  ~p  sin22Sin2Z 

.cic  a  b 


D  (l  +  a  )' 


9  4  u. 


^12 - ^Hr-P  sln?si“^ 

X12  D  (l  +  a2)2  a  b 


*  U42! 


q  a 
o 


~g  g  sin  —  sin  2Z 
D  (1  +  a2)2  a  b 


rf4  U  -  1  r»4  II  .  %  a 

V  U4l6 '2  7  ^62  ' 


D  (1  +  a2)2 


cos  22  cos  £Z 
a  d 


(8.13) 


4  14 

V  %> 6  "  2  V  U46i 


10a“ 


D  (1  +  a2)2 


cos  22  cos  «Z 
a  b 


and 


4  4 

V  wnll  =  V  wnl 1 


w. 


*nl£ 


- sin  22  sin  2Z 

D  (1  +  a2)2  a  b 
4 


4  * 

V  vn}£  =  tj  (l  + 


tq  a 

- - — w  'n  cos  22  cos  2Z 

3(1  +  a2)2  a  b 

4  «3oa2 

t7  v  *  _ _ -2__ 

V  n12  nfi+« 


— ~  22  sin  |2  Bin  2Z 

D  (1  +  a2)2  a  b 


w 


4 

7  w, 


k  (l  -yj  q  a2 

W  v  vrf-  D(1,ay  "“P-P 

4  ,  *%a3 


4q_  or 

- -g  g  g  cos  22  cos  2Z 

D  (l  +  a2)2  a  b 

r  ss  m  v  Mr 

n22  V  n22 


4 

W  “ 


*0“ 


(8.14) 
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The  remaining  equations  are  all  of  the  fom  ^^(l^w)  =  °>  with  their  solutions 
identically  equal  to  zero. 


The  above  equations  are  of  two  types : 


(1) 

and  (2) 


(  )  «  kl  sin  S  sin  p 


(3.15) 


V"  (  )  =  k2  cos  p  cos  p 


The  solutions  of  these  two  types  of  equations  are  (see  Reference  (25),  for  example): 


(1) 


ki a 


*  (1  +  ac) 


kx  «y 

2T2  sin  r  sin  b 


(8.1b) 


and  (2)  (  ) 


where 


n4(l  +  a2)2  .  A 


c  sin  Sin  WZ 
ran  a  b 


m»2, 4  n=2, 4 


lttan 


(» -  l)(n  -  1) 

Therefore,  the  complete  second  approximation  is : 


w  =  -T— 
it  D( 


1  +  (nn  +  nll} 


(1  +  or ) 


1  ...  _  ^  a 

2x2  n12  /,  .  _.2n2 


(1  +  ac)‘ 


(1  -20  a 

( 


(8.17) 


•  -"v)  a  «4  1 

T77F  (°66  * ns?)  +  (”22  +  “2|>  (TT^J 

7^^  [a<»i6  -  “jg> *  “3(”2s  -  ”^>1  I 


sin  **  sin  & 
a  b 


n*«2, 4  m=2, 4 


c  sin  —  sin^ 
mn  a  o 


84 


u  = 


n4D( 


[ 


7777  f  [<*i7  V  <  V  *2i> J  sl»  ¥ 

} 


^16+  2*62)a  +  (4o+  ^5i)q3  ] 


«0 

r-1 

/ 


m=2, 4  n=2, 4 


c  sin*^  sin”^ 
an  a  b 


£Z 

b 

(3.13) 


RECTAHGULAR  PLATE  WITH  TWO  OPPOSITE  EDGES  SIMPLY  SUPPORTED 
ADD  THE  OTHER  TWO  EDGES  CLAMPED 

As  in  the  previous  case,  a  sinusoidal  distribution  of  load  vill  be  assimed. 
For  the  coordinate  system  as  shown. 


y 

The  equation  for  the  load  is : 


q  =  qQ  sin  ^  cos  ^ 


(3.19) 


To  simplify  the  presentation,  the  arrangement  of  the  layers  will  be  assumed 
to  be  symmetrical  with  respect  to  the  middle  surface  of  the  plate.  An  extension 
to  a  general  arrangement  of  the  layers  can  be  readily  obtained  from  the  solution 
given. 


The  solution  is  assumed  to  consist  of  a  particular  and  a  homogeneous  part; 
i.e.,  w0  =  wi  +  w2.  The  particular  solution  is  chosen  so  as  to  satisfy  the  load¬ 
ing  condition,  whereas  the  homogeneous  solution  is  such  that  the  sum  of  the  two 
parts  satisfies  the  boundary  conditions.  The  particular  solution  is  taken  as 
that  for  the  case  when  all  edges  are  simply  supported: 


4 

*0  a 

D  (1  +  a2)2 


sm 


JtX 


cos 


2Z 

b 


(8.20) 
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where  a  =  a/b,  the  aspect  ratio. 

For  symmetry  with  respect  to  the  x  axis,  the  homogeneous  solution  can  be 
shown  to  be  (see,  for  example.  Reference  (25)  ): 


w  =  [a  Cosh  2Z  +  B  2Z  Sinh  22  1  sin  22 

2  [  a  a  a  J  a 


(8.21) 


where 


kSitihk 


TV  (l  +  a2)2  Sinh  ^  Cos  ^  ^ 


Cosh  §* 


TV  (1  +  a2)2  Sinh  ^  Cosh  %  +  ^ 


Therefore,  the  solution  for  the  first  approximation  is : 


w  =  Q  sin  —  [cos  22L  k  (2-  Sinh  Cosh  2Z  -  SZ  Cosh  2-  Sinh  £Z]  (8.22) 

o  a  I  b  23  33  a  a  d3  a  I 


where 


D  (l  +  a)‘ 


Sinh  §  Cosh  §_  +  §5 


Because  of  the  simplifying  assumption  regarding  symmetry  with  respect  to  the 
middle  surface,  the  only  differential  equations  for  the  corrective  teims  which 
need  to  be  considered  are  Eq  (4.l8).  Substitution  of  Eq  (8.22)  into  Eq  (4. l8) 
leads  to  the  following  equations  when  the  functions  are  expanded  in  suitable 
trigonometric  series : 

i)  r  ®  1 

4  ~  n  .  jtx  ny  ^  \  (nil)  n«y  I 

V  wnll  ’  51  T  sin  r  COS  b  ’  K  L  Cn  cos  b 

L  n-1,3  1 


where 


c'”11)  .(-1) 
n 


_ 8n  a 

x  (1  +  n2 
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where 


l  Jg.4  < -a  + °a>2 
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*m  +  n  at  ;  n.g>  4 

Sinh  ^  Cosh  J-  +  % 
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(-D2 


(nl5) 
na6n 

(m2  +  n2  a2)2  , 


(n66) 


v  <C£  =  (1 


where 
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Sinh  §j  Cosh  J-  +  % 
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■5 - 5—5  cos 

d.  <d.  \d  D 


where 


It  should  be  recalled  that  the  second  approximation  is  obtained  by  summing 
these  separate  solutions  times  the  appropriate  parameters : 

w  =  w  +  n„,w  +  n.^v  +  n^^w  ^  +  n^w  +  n_0w 

o  11  nil  16  nl6  66  n66  26  n26  22  n22 
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9-  OPTIMIZATION  OF  LAMINATE  CONFIGURATION 

The  question  of  obtaining  the  optimum  orientation  of  the  layers  in  a  structural 
laminate  can  be  considered  from  several  vievpoints  There  is  the  question  of 
whether  to  optimize  for  minimum  deformation  (or  minimum  volume  change  in  a  pressure 
vessel),  for  equal  stress  in  two  directions,  for  maximum  strength  (which  involves 
the  question  of  a  valid  theory  of  failure),  or  for  some  other  stress  or  deformation 
condition.  In  some  cases  two  conditions  may  be  satisfied  by  the  same  orientation, 
but  it  is  not  evident  that  this  should  be  true  in  general  Then,  the  particular 
application  for  which  the  optimum  structure  is  to  be  determined  may  be  classified 
as  to  whether  the  stresses  are  statically  determinate,  as  are  the  membrane  stresses 
in  a  thin-walled  pressure  vessel,  or  whether  the  stresses  are  statically  indeter¬ 
minate,  as  in  the  bending  of  a  plate.  Thus,  a  general  discussion  of  optimization 
would  become  quite  involved.  Only  one  aspect  of  the  problem  will  be  considered 
here. 


The  question  of  optimization  for  minimum  deflection  involves  only  a  study 
of  solutions  to  the  differential  equations  presented  in  the  first  part  of  this 
report.  This  discussion  of  optimization  will  be  limited  to  optimizing  for 
minimum  deflection  in  plate  problems .  Two  classes  of  problem  will  be  considered: 
the  cylindrical  bending  and  extension  of  long  rectangular  plates,  and  the  bending 
of  rectangular  plates  simply  supported  on  two  opposite  edges. 

CYLINDRICAL  BENDING  OF  LONG  RECTANGULAR  PLATES 

As  shovn  in  Section  5,  in  the  cylindrical  bending  of  long  rectangular  plates, 
the  lateral  deflection  may  be  expressed  in  the  following  form: 


w(x) 


"11 


(D11)£ 


F(x) 


(9-1) 


where  F(x)  may  be  obtained  from  the  corresponding  problem  for  the  deflection  of  a 
beam-  Similarly,  the  extension  of  a  long  rectangular  plate  may  be  expressed  in 
the  following  form: 


u  (x)  = 


A11  " 


Dn 


0(x) 


(9-2) 


It  can  be  seen  from  the  above  expressions  that  in  order  to  minimize  the 
deflection  or  extension  in  the  plate,  it  is  necessary  to  maximize  the  expression 
(D  *)2  (D  *)2 

D11 - A^i —  or  ^11 - —  »  respectively.  For  orientations  which  are 

symmetrical  with  respect  to  the  middle  surface  of  the  plate,  the  quantity  D1  *  is 
zero.  Therefore,  it  is  obvious  that  the  optimum  orientation  will  be  one  which 
is  symmetrical.  For  symmetrical  orientations,  the  deflection  or  extension  is 
minimized  by  maximizing  the  quantity  Di  or  A-|  i ,  respectively. 
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Fran  the  definitions  of  the  quantities  D-,  n  and  At  t  : 


Bn  -  f  L  bA*)  <>£  - 


k=l 


Au  ■  I  siik)  <\  -  Vl>' 


k=l 


(9-3) 


it  can  be  seen  that  in  order  to  maximize  these  quantities,  the  quantity  Cqq  must 
be  maximized  for  each  layer.  From  the  definition  of  C-y_: 


C±1  -  Cu  cos  0  +  2(0^  +lCgg)  sin 
the  conditions  for  maximizing  can  be 


sin  0=0 

y 

when 

C11> 

C12 

cos  0  =  0 

> 

when 

C22  * 

C11 

2  0 

(°12 

+  “«> 

C11 

tan  y  c  "■ 

"  C22 

> 

'  0  cos  0  +  C22  sin  0  (9.4) 

found  to  be: 


2C66 

“a  cu>  C22 

and 

C22>  °12  * 

(9-5) 

when 

C12*X66>  Cu 

and 

C12  +  “66  >  C22 

The  optimum  orientation  depends  on  the  relative  magnitudes  of  C]j  and  Cpp.  For 
Cxi  =  C22  >  the  optimum  orientation  is  45°,  as  vould  be  expected. 

RECTANGULAR  PLATES  SIMPLY-SUPPORTED  ON  TWO  OPPOSITE  EDGES 

From  the  solution  given  in  Section  8,  the  maximun  deflection  of  a  sinusoidally 
loaded  rectangular  plate  simply  supported  on  all  four  edges  can  be  shown  to  be: 


qo  a 


~  7  D(i  +  a2)2 


[  1  +  <»u  + 
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+  n12 
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(1  +a2)2 


+  (n, 


'66  +  n6# 


(1  -3/D)  a* 
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+  ^22  + 


°2l>  <l“a2)2] 


(9.6) 


where  a  =  a/b  is  the  aspect  ratio,  D  is  the  larger  of  Du  or  D12  +  2Dg£  or  D22> 
and  the  n.^  are  small  positive  quantities  defined  by: 
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nll  =  1 


n22  =  1 
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(9-7) 


From  this  solution,  it  can  be  seen  that  to  minimize  the  deflection  the 
quantity  D  should  be  maximized  and  the  ny  and  the  ny  should  be  minimized.  These 
two  effects  are  not  entirely  independent  since  a  change  in  orientation  which  de¬ 
creases  the  ny  will  also  decrease  D.  However,  since#for  arrangements  which  are 
symmetrical  with  respect  to  the  middle  surface  the  ny  will  be  zero,  the  optimum 
arrangement  will  be  one  which  is  symmetrical  with  respect  to  the  middle  surface. 

For  symmetrical  arrangements  of  the  layers,  Eq  (9*6)  can  be  written  as 
follows  (assuming  that  Djj_  >  Dp?) : 


when  D, 


>  +  D 


:i  +  a2)2[ 


Dn-(D^+2D66) 


Dll"D22 


(1+  a2)2 


a4 

(1+aV 


(9-3) 


when  (D^  +  2D66)  >  D^,  D  -  (D^  +  2D66) 

_  >  r 


*  (°i2+  ^  (1+ 


"a2)2  ; 


D-D,. 


(1+  a2)2 


(1+  a2)2 


(9.9) 


For  long  plates  (a  ■  0),  it  can  be  seen  that  the  first  of  these  equations  reduces 
to  one  similar  to  that  given  in  the  previous  discussion  of  long  plates,  whereas 
the  second  equation  is  different  in  form  since  two  terms  will  be  retained.  From 
the  previous  discussion,  it  should  be  realized  that  the  maximization  of  the  Dy 
depends  on  the  relative  magnitude  of  Cy  and  +  2C55  (Cy  can  always  be  defined 
to  be  greater  than  Cpp) ■ 

Rather  than  attempt  to  continue  in  general  terns,  a  series  of  numerical 
examples  will  be  studied.  First,  consider  the  material  in  the  layers  to  be  such 
that  Cqi  *  1.0C,  C1P  +  2Cgg  ■  C22  "  0.8C.  ( These  values  are  chosen  so  that  the 
ny  will  be  less  than  0.2,  and  their  squares,  products,  and  higher  powers  will  be 
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0  -  0°), 


C22  *  0,9C 

When  these  layers  are  oriented  so  that  Ci p  +  2Cgg  is  a  maximum  (i.e., 
it  can  be  shown  -that: 


=  0.3C 


c12  *  2C66  ’  1-°° 


Ca2  -  0.3C 


For  the  first  case  (0  **  *<-50),  Eq.  (9*8)  becomes: 


w  =  V! — r  i+ o.w — 

““  THT  (l  +  a2)2  L  (l  +  a2)2  J 

For  the  second  case  (0=0°),  Eq  (9-9)  becomes: 


(9.12) 


may  4 


D  ( 


1  f  .  .  0.2  0.2  a4  1 

1  +  a2)2  L  (1+a2)2  (l  +  a2)2J 


(9-13) 


As  an  indication  of  the  effect  of  the  aspect  ratio,  a  ,  the  maximum  deflection 
has  been  evaluated  from  Eqs  (9.10),  (9-H),  (9-12),  and  (9-13)  for  the  limiting 
conditions  a  *  0  (i.e.,  the  long  plate)  and  a  -  1  (i.e.,  the  square  plate).  The 
results  can  be  summarized  as  follows  : 
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(©  =  o°) 

1. l6 

(0  -  ^5°) 

0.287 

(0  -  0°) 
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(0  =  45°) 

=  0.3c 

C12  +  2C66  =  1,0(3 

C22  =  0.8C 

1.111 

(©  -  ^5°) 

1. 20 

(0  =  0°) 

0.309 

(0  =  ^5°) 

0.275 

(©  -  o°) 

From  the  left  hand  side  of  this  table,  it  can  be  seen  that  for  long  plates 
(a  =  0)  the  layers  should  be  oriented  so  as  to  maximize  D-, ^  (the  same  conclusion 
indicated  previously).  From  the  right  hand  side  of  this  table,  it  can  be  seen 
that  for  square  plates  (a  =  l)  the  optimum  orientation  tends  to  be  one  which 
maximizes  Dtp  +  2 D^g,  although  the  difference  between  the  values  in  the  first 
line,  0.237  and  O.276,  is  within  the  accuracy  to  be  expected  of  the  perturbation 
method  of  solution,  (using  orthotropic  plate  theory,  which  is  applicable  in  this 
case,  it  can  be  shown  that  the  value  O.287  should  be  0.278). 

Interpreting  these  results  in  terms  of  the  angle  of  orientation,  it  appears 
that  when  Cj_j_  is  greater  than  C-|  p  +  2Cg6  the  angle  of  orientation  of  each  layer 
should  be  zero.  On  the  other  hand,  when  C12  +  2C66  is  greater  than  C11,  it  appears 
that  for  long  plates  the  angle  of  orientation  should  be  45  degrees  (for  C-j^  =  Cpp), 
whereas  for  square  plates  the  angle  of  orientation  should  be  zero.  The  aspect 
ratio  for  which  the  optimum  angle  of  orientation  switches  from  45  degrees  to  0 
degrees  may  be  determined  by  setting  Eq  (9-8)  equal  to  Eq  ( 9 - 9 ) *  Such  a  computa¬ 
tion  for  the  material  constants  indicated  in  the  second  line  of  the  above  table 
gives  an  aspect  ratio  of  0.45.  However,  the  accuracy  of  such  a  computation  is 
somewhat  questionable  because  of  the  approximate  nature  of  the  method  of  solution. 
The  results  may  be  checked  by  more  accurate  methods  of  solution  (a  refined 
perturbation  method,  orthotropic  plate  theory,  or  anisotropic  plate  theory), 
since  only  special  orientations,  for  which  these  methods  are  applicable,  need  to 
be  considered. 

For  the  case  of  a  rectangular  plate  simply  supported  on  two  opposite  edges 
and  clamped  on  the  other  two  edges,  a  study  of  optimum  orientations  similar  to 
that  given  above  could  be  made.  The  conclusion  that  the  optimum  arrangement  of 
the  layers  is  one  which  is  symmetrical  with  respect  to  the  middle  surface  of  the 
plate  would  be  found  to  apply  in  this  case  also.  Furthermore,  the  solution  can 
be  shown  to  reduce  to  be  similar  in  form  to  that  for  long  rectangular  plates  when 
a  «  0.  However,  it  can  be  seen  from  the  solution  given  in  Section  8  that  it  will 
require  extensive  numerical  calculations  in  order  to  arrive  at  specific  conclusions 
regarding  the  optimum  orientations  for  particular  aspect  ratios. 

EFFECT  OF  THE  NUMBER  OF  LAYERS  FOR  A  GIVEN  TOTAL  THICKNESS 

As  indicated  in  the  previous  discussion,  for  many  cases  the  deflection  is 
minimized  by  orienting  all  of  the  layers  in  one  direction.  In  these  cases,  if 
the  effect  of  the  glue  layers  is  neglected  and  if  the  layers  are  all  of  the  same 
material,  there  will  be  no  effect  of  the  number  of  layers.  If  the  effect  of  the 
glue  layers  is  considered,  there  will  be  a  slight  increase  in  the  deflection  as 
the  number  of  layers  increases. 
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.ART  III:  DETERMINATION  OF  MATERIAL  PROPERTIES,  FRACTURE  CRITERIA 
INTRODUCTION 

This  part  of  the  report  contains  a  section  describing  a  possible  procedure 
which  can  be  followed  to  determine  the  elastic  properties  of  an  individual 
orthotropic  lamina.  Such  properties  are  necessary  for  use  in  the  plate  and  shell 
theory  developed  in  Part  I.  The  final  section  presents  a  brief  description  and 
the  results  of  a  limited  number  of  uniaxial  crack  propagation  tests  of  laminated 
foil  type  laminates. 
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10.  DETERMINATION  OP  ELASTIC  COEFFICIENTS  OF  INDIVIDUAL  LAMINAS 


As  described  in  Part  I,  Section  1,  four  independent  elastic  coefficients 
are  required  to  specify  the  stress-strain  relations  for  an  orthotropic  lamina. 

These  are  conveniently  taken  to  be  tensile  moduli  of  elasticity  E^  and  Ep,  shear 
modulus  of  elasticity  Glt  and  either  value  of  Poisson'6  Ratio  VpT  or 
Since  it  appears  to  be  difficult  to  measure  the  shear  modulus  of  foil  directly, 
em  indirect  method  must  be  employed.  One  possible  approach  is  to  perform  tensile 
stress -strain  tests  of  a  foil  loaded  uniaxially  first  in  the  L  direction,  then 
the  T  direction,  and  finally  in  one  other  direction  (say  45  degrees).  If  the  first 
two  tests  El,  Ej  and  3^lt  or  ^TL  can  be  measured,  the  value  of  GLT  can  be  calcu¬ 
lated  from  the  third  test.  Ihe  details  of  this  procedure  follow. 


Following  the  steps  outlined  in  Part  I,  Section  1,  it  is  possible  to  write 
strain-stress  equations  for  an  orthotropic  lamina  referred  to  arbitrary  axes 
x,  y.  Eq  (1.2)  then  takes  the  more  general  fom 
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The  relationships  between  the  Spj  and  Spj  and  the  angle  ©  are  given  by 

cos4  ©  +  (2S10  +  S ,A  sin2  ©  cos2  ©  +  S0_  sin4  © 

11  11  '  12  66'  22 

(cos4  ©  +  sin4  ©)  +  (S1;L  +  S22  -  S66)  sin2  ©  cos2  ©  (10.2) 

Spp  ■  Sg2  cos4  ©  +  (2S^2  +  Sgg)  sin2  ©  cos2  ©  +  sin4  © 

Sgg  =  Sgg  (cos4  ©  +  sin4  ©)  +  2(2S1:l  +  2S22  -  -  Sgg)  sin2  ©  cos2  © 


Sl6  -  (2SX1  -  23^  -  Sgg)  cos3  ©  sin  ©  -  (2S, 
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32g  -  (2Si;l  -  2S^2  -  Sgg)  cos  ©  sin3  ©  -  (2S22  -  2S1P  -  S^)  cos3  ©  sin  © 


12  66 ; 


The  constants  Spp  and  Sop  can  be  obtained  from  tension  tests  with  ©  ■  0C 
and  90°  respectively,  since  Spp  becomes : 

for  ©  «*  0° 
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(10.3) 
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and  for  ©  =  90° 
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22 


©  =  45° 


The  remaining  constants  to  be  evaluated  are  Sip  and  S55;  however,  the  only  value 
which  can  be  obtained  directly_by  tension  tests  is  the  value  (2S ^  +  Sgg).  If 
©  =  45°  in  the  expression  for  then 


<2S12* 


s66> 


4s 


11 
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I©  =  45° 


(10.5) 


To  obtain  Sgg,  Sjn  must  be  experimentally  measured.  This  involves  the 
measurement  of  Poissorrs  Ratio.  A  possible  procedure  would  involve  means uring 
V  LT  and  VTL,  from  which 


2S12  =  S12  +  S21 


*^LT  +  ~^TL  1 

'  et  el  J 


Substituting  Eq  (10.6)  and  Eq  (l.  3)  into  Eq  (10.5)  gives 


S66 
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E45° 


(10.6) 


(10.7) 


In  this  way  the  required  elastic  coefficients  may  be  determined. 

Preliminary  tensile  tests  on  1100  -  H19  aluminum  foil  of  0.0035  inch  thickness 
indicate  that  mechanical  property  values  can  be  obtained  as  described  above.  To 
obtain  complete  information,  one  must  conduct  tension  tests  in  the  longitudinal, 
transverse  and  45  degrees  directions  from  which  moduli  of  elasticity  E^,  ET, 
and  E45  degrees  8,8  well  Poisson's  Ratios  and  3/TL  can  be  determined. 

The  photogrid  method  which  has  been  successfully  applied  to  a  number  of  deformation 
problems  in  metal  sheets  can  be  applied  to  the  problem  of  determining  the  fore¬ 
going  parameters. 
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11.  CRACK  PROPAGATION  IN  ALUMINUM  FOIL  LAMINATES  * 

Crack  propagation  tests  of  4  x  10  inch  specimens  load  in  uniaxial  tension  vas 
performed  to  study  the  fracture  characteristics  of  foil-type  laminates.  The  material 
used  for  making  the  laminated  specimens  vas  1100  H-ig  aluminum  foil,  0.0037  inch 
thick.  The  ultimate  tensile  strength,  0.2  per  cent  tensile  yield  strength  and 
elastic  modulus  were  obtained  from  1x4  inch  specimens  tested  on  an  lustron 
(Model  TT-EL)  tensile  testing  machine.  The  results  of  the  tests  are  shown  in 
Table  4.  Loading  in  the  grain  direction,  across  it  as  veil  as  at  45  degrees  to 
the  grain  direction  of  the  foil,  vas  done  to  determine  the  anisotropy  of  the 
aluminum  foil. 

The  laminated  specimens  vere  made  from  4  x  10  inch  sheets  of  1100  H-lg 
aluminum  foil  bonded  vith  epoxy  resin  FM-47 • Measurements  of  the  glue-line 
thicknesses  shoved  that  they  vere  between  0.00004  to  0.00016  inch;  these  glue  lines 
vere  achieved  by  applying  a  3:2  thinner  to  epoxy  resin  mixture  vith  a  special  spray 
gun.  A  typical  specimen  configuration  is  shown  in  Figure  1,  and  as  can  be  seen, 
the  initial  crack  vas  made  by  drilling  a  1/6- inch  hole  and  saving  a  0.00S  inch 
slot  of  desired  length  through  the  entire  thickness  of  the  laminated  specimens. 

Two,  four,  six,  and  ten-ply  laminated  sheets  vith  the  foil  grain  direction  both 
parallel  and  perpendicular  to  the  direction  of  the  applied  tensile  load  vas 
tested.  The  specimens  having  initial  crack  lengths  xQ  of  0.5,  0.6,  0.6,  and  1.0 
inches  vere  tested  in  a  Tate-Baldvin  60,000  pound  tensile  testing  machine.  A 
scale  of  0.01  inch  least  count  vas  used  to  measure  crack  extension  during  the  slow 
load  application  until  unstable  crack  lengths,  i.e.  sudden  fracture,  occurred. 

The  experimental  results  are  summarized  in  Table  5  and  are  plotted  in  Figures 
2,  3,  4,  and  5  for  2-,  4-,  6-,  and  10-ply  laminated  specimens  respectively.  The 
plots  of  gross  area  stress  vs.  crack  length  show  the  crack  initiation  and  failure 
lines  for  specimens  of  varying  crack  length.  During  the  course  of  the  experiments, 
it  vas  observed  that  upon  reaching  maximum  load  a  brief  discontinuous  crack  growth 
took  place  without  the  usual  load  drop  experienced  in  tests  of  monolithic  specimens. 
Although  brief  in  nature,  this  critical  propagation  characteristic  illustrates  the 
temporary  ability  of  the  laminate  to  arrest  propagation  until  loads  beyond  the 
maximum  critical  load  are  applied.  In  monolithic  specimens,  the  crack  extension 
after  critical  conditions  are  reached  at  loads  below  the  maximum  load.  Both  the 
monolithic  sheets  and  foil  laminates  failed  after  formation  of  deep-necked  bands 
in  the  plane  of  the  sheet.  The  bands  in  the  monolithic  sheet  specimens  formed 
at  45  degrees  to  the  applied  tensile  load,  the  bands  in  the  laminate  specimens 
vere  not  as  deep  and  formed  on  a  fracture  plane  normal  to  the  applied  load.  The 
lesser  necking  of  the  laminates  in  the  fracture  region  means  that  reduction  of 
cross  sectional  area  is  less  than  in  monolithic  sheets  and  consequently  the  load 
drop  is  also  reduced.  The  values  of  net  area  initiation  and  fracture  stress 
shown  in  Table  4  show  that  crack  propagation  initiates  at  a  stress  far  below 
(less  than  60  per  cent)  the  ultimate  strength  of  the  material.  The  net  area 
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stress  approaches  the  ultimate  strength  when  the  critical  crack  length  is  reached. 

Figures  6a  and  6b  shov  plots  of  final  crack  length  as  a  function  of  the 
number  of  foils  in  the  laminate.  The  ability  of  10-ply  specimens  to  withstand 
longer  final  crack  length,  for  particular  initial  crack  length,  them  2-,  4-,  and 
6-ply  laminates  is  illustrated  and  can  be  directly  correlated  to  the  10-ply 
laminate's  greater  load  carrying  capacity  as  shown  in  Table  4.  However,  direct 
comparison  between  2-  and  10-ply  specimens  should  be  done  only  with  consideration 
of  the  more  pronounced  buckling,  i.e.,  lower  fracture  strength,  of  the  2-ply 
specimens.  The  linear  relationship  between  final  and  initial  crack  length  for 
different  laminates  is  shown  in  Figures  7a  and  Tb  for  with-grain  and  cross-grain 
loading  conditions. 

Gross  area  initial  stress  and  maximum  gross  area  stress  as  functions  of 
initial  crack  length  are  shown  in  Figures  3  and  9>  respectively.  These  data  have 
been  replotted  on  log- log  coordinates  in  Figures  10  and  11  to  test  the  theoretical 
relationship  that  the  maximum  gross  area  stress  will  vary  inversely  with  the  square 
root  of  the  initial  crack  length.  Figures  11a  and  lib  show  that  curves  of  slope 
0.5  can  be  fitted  to  the  experimental  data  of  each  specimen  group  and  the  power 
relationships  are  tabulated  in  Table  5*  However,  a  similar  attempt  to  relate  the 
gross  area  initiation  stress  to  initial  cracklength,  shown  in  Figures  10a  and  10b, 
gives  no  general  relationship  as  indicated  in  Table  5>  where  the  exponents  vary 
from  0.14-7  to  O.663.  The  gross  area  crack  initiation  stresses  and  the  maximum 
gross  area  stresses  as  functions  of  final  crack  length  xf  before  sudden  fracture 
are  plotted  in  Figures  12  and  13-  Since  the  final  crack  length  is  related  to  the 
initial  one  as  shown  in  Figure  7>  the  resultant  decrease  in  strength  with  final 
crack  length  as  6hown  in  Figures  12  and  13  follows  the  patterns  of  the  similar 
plots  with  respect  to  the  initial  crack  length. 

Of  particular  interest  is  the  Griffith- Irwin  (26)  fracture  criterion  dtf/dA, 
the  dissipation  rate  of  plastic  work  during  fracture,  shown  as  a  function  of  initial 
crack  length  in  Figures  l4a  and  14b.  It  can  be  observed  that  the  values  of  dW/ dA 
remain  relatively  constant  for  each  laminated  specimen  group.  The  6-  and  10 -ply 
specimens  show  a  substantially  greater  value  of  dW/dA  than  do  the  2-  and  4-ply 
specimens.  Since  higher  dissipation  rates  are  associated  with  slower  crack  pro¬ 
pagation  rates  before  sudden  fracture,  it  may  be  concluded  that  additional  foil 
layers  would  be  beneficial.  However,  as  shown  in  Figure  l4a  the  dissipation  rate 
would  be  increased  only  slightly.  The  calculated  values  of  dtf/dA  for  each  specimen 
and  the  average  values  for  each  specimen  group  are  shown  in  Table  5-  These 
calculated  values  incorporate  the  correction  for  the  ratio  of  initial  crack  length 
to  specimen  width.  It  is  noteworthy  that  the  average  value  of  dW/dA  for  any 
specimen  group  when  substituted  in  the  Griffith-Irvin  equation  for  maximum  gross 
area  stress  will  yield  essentially  the  same  constant  shown  in  Table  5  for  the 
empirical  power  function  between  maximum  gross  area  stress  and  initial  crack 
length. 

For  example,  the  empirical  equation,  as  shown  in  Table  5>  for  2-ply  laminates 
loaded  in  the  cross -grain  direction  is 


0 

max 


“  *3-°  xq 


-0.475 


(ksi) 
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The  average  value  of  d///dA  =  30.6  together  vith  an  elastic  modulus  of  10?  psi 
when  substituted  in  the  Griffith- Irvin  equation  gives 


a  =  13.9  x 
max  o 


-0.500 


(ksi) 


vhich  is  sufficiently  close  to  pemit  failure  predictions  based  on  dissipation 
rate  values. 


Table  4 


Material  Properties  of  1100  E-19  Aluminum  Foil  Sheets 


Ultimate  Tensile 
Strength 
ksi 


0.2$  Tensile 
Yield  Strength 
ksi 


Modulus  of 
Elasticity 

psi  x  10-6 


Tensile  Load  with 

Grain  Direction  of  Foil  31-10  28.69  9-87 

(with  grain) 


Tensile  Load  at  ^5°  27  38 

to  Grain  Direction  of  Foil 


9-15 


Tensile  Load  at  90° 

to  Grain  Direction  of  Foil  23.50  2b. 10  9-33 

(cross  grain) 
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Table  5 

Summary  of  Results  of  Crack  Propagation  Tests 
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10 PLY  LAMINATE  0.500  0.895  0.1340  20.45  15.23  52.7  -o.«i8  _  _  ■ 

LOADING:  WG  1.030  1  500  0.1388  16.15  12.98  51.9  J  xi 


Foil  Laminate  Specimen 


Figure  2a.  Gross  Area  Stress  vs  Crack  Length  of  2-Ply  1100  H-19  Aluminum  Foil  Laminates 
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CRACK  LENGTH,  X (INCHES) 


Figure  2b.  Gross  Area  Stress  vs  Crack  Length  of  2-Ply  1100  H-19  Aluminum  Foil  Laminates 
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CRACK  LENGTH,  x  (INCHES) 
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CRACK  LENGTH  ,  x  (INCHES) 


Figure  3b.  Gross  Area  Stress  vs  Crack  Length  of  4-Ply  1100  H-19  Aluminum  Foil  Laminates 
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CRACK  LENGTH,  X  (INCHES) 
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Figure  4b.  Gross  Area  Stress  vs  Crack  Length  of  6-Ply  1100  H-19  Aluminum  Foil  Laminates 
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CRACK  LENGTH  ,  x  (INCHES) 


Figure  5a.  Gross  Area  Stress  vs  Crack  Length  of  10-Ply  1100  H-19  Aluminum  Foil  Laminates 


112 


ENGTH,  x  (INCHES) 


Figure  5b.  Gross  Area  Stress  vs  Crack  Length  of  10-Ply  1100  H-19  Aluminum  Foil  Laminates 


CRACK  LENGTH  ,  x  (INCHES) 


NUMBER  OF  FOILS  IN  LAMINATE 

Figure  6a.  Final  Crack  Length  vs  Number  of  Ply  in  1100  H-19  Aluminum 
Foil  Laminates 
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FINAL  CRACK  LENGTH,  xf  (INCHES) 


I00H-I9  ALUMINUM  FOIL  LAMINATES 
LOADING:  CROSS  GRAIN 


123456789  10 

NUMBER  OF  FOILS  IN  LAMINATE 

Figure  6b.  Final  Crack  Length  vs  Number  of  Ply  in  1100  H-19 
Aluminum  Foil  Laminates 
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FINAL  CRACK  LENGTH,  xf  (INCHES) 


INITIAL  CRACK  LENGTH,  x0(lNCHES) 

Figure  7a.  Final  Crack  Length  vs  Initial  Crack  Length  in 
1100  H-19  Aluminum  Foil  Laminates 
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Figure  7b.  Final  Cr 
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Figure  8a.  Gross  Area  Initiation  Stress  vs  Initial  Crack  Length 
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ROSS  AREA  INITIA 


I 


| 


INITIAL  CRACK  LENGTH,  x0  ( INCHES) 

Figure  8b.  Gross  Area  Initiation  Stress  vs  Initial  Crack  Length 
in  1100  H-19  Aluminum  Foil  Laminates 
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MAXIMUM  GROSS  AREA  STRESS, 


INITIAL  CRACK  LENGTH,  x0  ( INCHES) 

Figure  9a.  Maximum  Gross  Area  Stress  vs  Initial  Crack  Length 
in  1100  H-19  Aluminum  Foil  Laminates 
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turn  Gross  Area  Stress  vs  Initial  Crack  Length 
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GROSS  AREA  INITIATION  STRESS,  (77  (KS 1 ) 
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INITIAL  CRACK  LENGTH,  x0  (INCHES) 


Figure  10a.  Log-log  Plot  of  Gross  Area  Initiation  Stress  vs  Initial 
Crack  Length  in  1100  H-19  Aluminum  Foil  Laminates 
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GROSS  AREA  INITIATION  STRESS,  Oj  IKSI) 


0-2  0.4  0.6  0.8  1.0  2.0  4.0  6.0 

INITIAL  CRACK  LENGTH.  x0  (INCHES) 


Figure  10b.  Log-log  Plot  of  Gross  Area  Initiation  Stress  vs  Initial 
Crack  Length  in  1100  H-19  Aluminum  Foil  Laminates 
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MAXIMUM  GROSS  AREA  STRESS,^  (KSI ) 
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INITIAL  CRACK  LENGTH,  x0  (INCHES) 


Figure  11a.  Log-log  Plot  of  Max.  Gross  Area  Stress  vs  Initial 

Crack  Length  in  1100  H-19  Aluminum  Foil  Laminates 
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MAXIMUM  GROSS  AREA  STRESS.  0^AX  (KSI) 
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INITIAL  CRACK  LENGTH,  x0  (INCHES) 


Figure  lib.  Log-log  Plot  of  Max.  Gross  Area  Stress  vs  Initial 

Crack  Length  in  1100  H-19  Aluminum  Foil  Laminates 
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GROSS  AREA  INITIATION  STRESS,  <rT  (KSI ) 


FINAL  CRACK  LENGTH,  xf  (INCHES) 

Figure  12a.  Gross  Area  Initiation  Stress  vs  Final  Crack 

Length  in  1100  H-19  Aluminum  Foil  Laminates 
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GROSS  AREA  INITIATION  STRESS,  <TX  (KSI) 


FINAL  CRACK  LENGTH,  xf  (INCHES) 

Figure  12b.  Gross  Area  Initiation  Stress  vs  Final  Crack 

Length  in  1100  H-19  Aluminum  Foil  Laminates 
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MAXIMUM  GROSS  AREA  STRESS 


Figure  13a.  Max.  Gross  Area  Stress  vs  Final  Crack  Length 
in  1100  H-19  Aluminum  Foil  Laminates 
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MAXIMUM  GROSS  AREA  STRESS,  <rmQX  ( KS I ) 


FINAL  CRACK  LENGTH,  xf  ( INCHES) 


Figure  13b.  Max.  Gross  Area  Stress  vs  Final  Crack  Length 
in  1100  H-19  Aluminum  Foil  Laminates 
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INITIAL  CRACK  LENGTH,  x0( INCHES ) 

Figure  14a.  Plastic  Work  Dissipation  Rate  vs  Initial  Crack 
Length  in  1100  H-19  Aluminum  Foil  Laminates 
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INITIAL  CRACK  LENGTH,  x0(  INCHES) 

Figure  14b.  Plastic  Work  Dissipation  Rate  vs  Initial  Crack 
Length  in  1100  H-19  Aluminum  Foil  Laminates 
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